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1.0  Executive  Summary 

During  the  final  stages  of  this  project  we  encountered  some  delays  in  securing  a  no- 
cost  extension  that  was  requested  in  order  to  complete  additional  analyses  that  were  not 
originally  anticipated.  These  analyses  include  one  additional  radiometric  age  and 
geochemical  analyses  of  lava  flow  units.  This  Final  Technical  Report  serves  to  convey 
the  results  of  the  efforts  of  Tasks  1, 2, 3  and  4  (Table  1.1).  Because  Task  4  represents 
Integration  and  Interpretation,  the  results  contained  herein  largely  reflect  efforts  imder 
this  task.  It  is  expected  that  two  manuscripts  for  publication  will  stem  directly  from  the 
efforts  of  this  project.  An  additional  4  manuscripts  for  publication  are  expected  to  result 
from  efforts  that  set  the  stage  for  the  work  conducted  under  fius  contract. 

Seismogenic  strain  at  the  Wild  Horse  Mesa  suggests  that  deformation  is  partitioned 
with  depth,  whereas  the  time-integrated  record  of  shearing  suggests  a  much  finer 
spatial  scale  of  deformation  partitioning.  Active  deformation  is  characterized  by 
subhorizontal  east-west  maximiim  stretching  both  above  and  below  a  deptii  of  5  km. 

At  the  shallower  crustal  levels  maximum  shortening  is  subvertical  and  associated  wifii 
crustal  thinning,  whereas  at  die  deeper  crustal  levels  maximum  shortening  is 
subhorizontal  and  associated  with  horizontal  dextral  shearing.  Both  levels  record  plane 
strain  and  evidence  for  block  rotation  at  rates  different  than  for  the  large-scale 
continuum.  The  time-integrated  record  of  shearing  at  upper  Wild  Horse  Mesa  indicates 
that  horizontal  dextral  shearing  with  subhorizontal  west  southwest-east  northeast 
maximum  stretching  is  the  dominant  strain  geometry  recorded.  Lesser  crustal 
tiiickening  and  crustal  thinning  are  recorded  as  well,  and  the  faults  that  accommodate 
these  strains  are  spatially  intermixed  with  the  more  dominant  strike-slip  faults  that 
accommodate  horizontal  shearing.  Statistical  results  at  90%  confidence  suggest  that  the 
horizontal  shearing  is  characterized  by  prolate  strain  and  non-zero  relative  vorticity. 

The  former  is  consistent  with  transtension  and  the  latter  is  in  accord  with  clockwise 
vertical-axis  rotation  (viewed  down;  here  and  throughout  this  document  vertical-axis 
rotation  sense  is  given  looking  down)  of  crustal  blocks. 

Pliocene  lavas  and  sediments  of  the  Wild  Horse  Mesa  exhibit  this  clockwise  vertical- 
axis  rotation  of  fault  boimded  blocks.  We  measure  these  finite  rotations 
paleomagnetically,  relative  to  two  different  reference  frames.  At  two  localities  we 
average  secular  variation  through  sedimentary  sections  to  reveal  rotation  or  its  absence 
relative  to  paleogeographic  north.  Where  sediments  are  lacking  we  sampled 
geochemically-identifiable,  areally-extensive  lavas  to  measure  the  rotation  of  localities 
relative  to  one  another.  In  the  Coso  Wash  the  fanglomerate  member  of  the  Coso 
formation  exhibits  approximately  21.6®±10.7°  clockwise  vertical-axis  rotation  since  ca.  3 
Ma.  In  the  Wild  Horse  Mesa,  3.1-3.4  Ma  lavas  exhibit  about  12.0°±4.6°  (weighted  mean) 
clockwise  vertical-axis  rotation  relative  to  outcrops  near  Dead  End  Cabin.  In  aggregate 
the  data  suggest  that  the  Wild  Horse  Mesa  and  Coso  Wash  hosts  a  segment  of  3\e 
Eastern  California  shear  zone  which  accomodates  dextral  shear  by  vertical-axis  rotation 
of  blocks  above  a  deep-seated  right-lateral  fault.  Using  the  measured  fiiute  rotation  and 
a  piimed  block  model,  we  calculate  a  minimum  of  1.5  km  of  dextral  offset  across  the 
Wild  Horse  Mesa  oriented  N30°W  since  ca.  3.1  Ma. 

At  adjoining  sites  in  the  greater  Coso  Geothermal  Area  and  the  Darwin  Plateau 
other  workers  suggest  no  vertical-axis  rotation.  However,  the  former  exhibits  down  to 
the  west  tilting.  Thus,  style  of  block  rotation  appears  partitioned:  vertical-axis  rotation 
in  the  Wild  Horse  Mesa  and  horizontal-axis  rotation  (tilting)  in  the  geothermal  area. 
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suggesting  partitioning  of  the  master  faults  controlling  the  near-surface  rotation  and 
tilting. 


2.0  Introduction 

2.1  Project  Background 

The  overarching  goal  of  this  study  is  to  establish  an  understanding  of  the  kinematic 
history  of  the  Wild  Horse  Mesa  and  adjacent  crust  in  the  vicinity  of  the  Coso 
geothermal  field.  The  approach  adopted  takes  advantage  of  the  fact  that  die  basement 
of  the  Wild  Horse  Mesa  is  yoimg  (late  Cenozoic),  faulted,  well-exposed  and  actively 
accommodating  seismogenic  strain.  In  particular,  we  have  employed  geochemistry, 
paleomagnetism  and  geochronology  to  constrain  volcano-stratigraphy,  the  ages  of  lava 
flows  and  the  finite  rotation  of  flows  after  cooling.  Seismic  events  have  been  used  to 
characterize  active  strain  and  outcrop-scale  brittle  faults  have  been  used  to  examine  the 
time-integrated  record  of  shearing  deformation.  The  net  result  of  this  work  is  (1)  an 
imderstanding  of  the  geometries  and  ages  of  the  basalt  and  andesite  flows  that 
dominate  the  volcanic  basement  of  the  Wild  Horse  Mesa,  and  (2)  a  comparison  between 
the  instantaneous  non-recoverable  strain  reflected  in  earthquake  data  and  the  time- 
integrated  brittle  shearing  reflected  in  the  faults  that  deform,  and  produce  finite  block 
rotations  of  the  Wild  Horse  Mesa. 

2.2  Tectonic  Setting 

The  Coso  geothermal  field  is  approximately  centered  at  a  right-releasing  step  in  the 
diffuse  strike-slip  boundary  between  the  Sierran  microplate  and  the  southern  Basin  and 
Range  province  to  the  soulheast  (Unruh  et  al.  2002).  The  Wild  Horse  Mesa  lies 
immediately  east  of  the  geothermal  field  and  is  well  situated  to  compare  seismogenic 
(i.e.,  instantaneous)  strain  with  the  time-integrated  strain  recorded  by  outcrop-scale 
brittle  faults  and  deflection  of  paleomagnetic  decimations  within  the  late  Cenozoic  lava 
flows  that  compose  the  mesa.  This  location  is  the  site  of  faulting  that  probably  links 
horizontal  shearing  (i.e.,  strike-slip  dominated  strain)  in  the  Indian  Wells  Valley  with 
the  horizontal  shearing  north  of  Wild  Horse  Mesa.  The  latter  shearing  presumably 
links  with  the  Owen  Valley  fault  zone  in  the  area  of  the  1872  Owens  Valley  earthquake 
rupture. 

The  lava  flows  that  form  Wild  Horse  Mesa  are  transected  by  an  anastamosing 
network  of  normal  fault  scarps  (i.e.,  the  step-faulted  terrane  of  Duffield  et  al.  1980) 
indicating  a  history  of  crustal  thinning  within  a  strike-slip  fault  system.  The  dominant 
fault  scarps  are  northwest-facing  in  northern  Wild  Horse  Mesa  and  merge  into  east-  and 
west-facing  scarps  in  southern  Wild  Horse  Mesa.  The  lava  flows  cut  by  these  faxolts  are 
ca.  3.6  Ma  or  yoimger.  Stratigraphic  and  structural  observations  (see  Monastero  et  al. 
2002  for  a  review)  indicate  that  significant  northwest  dextral  shearing  initiated  in  this 
region  between  ca.  3.5  and  2  Ma,  and  geodetic  data  (Dixon  et  al.  2000)  indicate  tiiat 
approximately  22%  of  Pacific-North  America  plate  motion  is  presently  accommodated 
here  (~llmm/yr).  In  total,  the  Wild  Horse  Mesa  provides  a  unique  opportunity  to 
examine  die  kinematic  evolution  of  a  youthful,  transtensional  system. 

2.2.1  Seismicity 

Seismicity  in  the  vicinity  of  the  Coso  Geothermal  field  is  heterogeneously 
distributed  (focal  mechanisms  and  locations  for  the  seismic  events  used  in  this  study 
include  events  provided  by  E.  Hauksson).  Three  notable  clusters  of  seismic  events,  in 
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the  eastern  part  of  the  Indian  Wells  Valley,  beneath  the  southern  part  of  Wild  Horse 
Mesa  and  immediately  nortih  of  Wild  Horse  Mesa,  have  been  used  to  delineate  to  a  first- 
order  the  geometry  of  the  right-releasing  stepover  that  seperates  the  Sierran  microplate 
from  the  southern  Basin  and  Range  Province  (Figure  2.1).  The  geometry  of  this 
boimdary  will  no  doubt  be  refined  as  the  catalog  of  well-located  seismic  events  becomes 
larger. 

The  fact  that  the  Wild  Horse  Mesa  occupies  a  position  within  the  right-releasing 
stepover  in  the  locus  of  strike-slip  deformation  east  of  the  Sierra  Nevada  is  apparent  on 
the  basis  of  seismogenic  strain.  Seismicity  from  0-5  km  accommodates  crustal 
thinning  whereas  that  from  5-8  km  accommodates  horizontal  shearing  (Unruh  et  al. 
2002).  Additionally,  the  former  seismic  events  are  somewhat  spatially  disorganized, 
whereas  the  latter  define  a  planar  seismic  source  zone  that  strikes  subparallel  to  the 
trajectory  of  Sierran-North  America  motion  (Dixon  et  al.  2000;  Unruh  et  al.  2002). 

2.2.2  Volcano-stratigraphy 

The  Coso  volcanic  field  is  composed  of  multiple  generations  and  centers  of  late 
Cenozoic  alkalic  to  subalkalic  (Novak  and  Bacon,  1986)  volcanism  spanning  the  last  6 
Ma  that  Duffield,  Bacon  and  coworkers  have  studied  extensively  (e.g.  Bacon  et  al.,  1982; 
Duffield  and  Bacon,  1981).  We  retain  the  terminology  of  Duffield  et  al.  (1980),  defined 
by  weight  %  Si02:  basalt  48-54%,  andesite  54-60%,  and  dacite  60-65%. 

Volcanism  in  the  Coso  Range  began  with  basaltic  through  dacitic  lavas  erupted  5-6 
Ma  at  the  present  southeastern  edge  of  Owens  Lake.  The  bulk  of  the  rocks  studied  here 
derive  from  the  Wild  Horse  Mesa  (Figure  2.2)  which  is  composed  of  basaltic  through 
dacitic  lavas  previously  tiiought  to  have  erupted  between  3  and  4  Ma  into  an  oblique 
puU-apart  basin.  This  was  the  most  voluminous  phase  of  eruption  in  the  Coso  Range 
during  the  late  Cenozoic.  Following  tins,  dacite  flows  through  ryholite  pyroclastics 
erupted  between  2.5  and  3.5  Ma  at  flie  western  edge  of  the  range  near  Haiwee 
Reservoir.  Later,  a  handful  of  basalt  and  andesite  eruptions  occurred  around  the 
Tertiary-Quaternary  boundary  near  Renegade  Canyon  on  the  Wildhorse  Mesa.  Since 
1.1  Ma,  bimodal  basalt  and  small-volume  rhyolite  volcanism  has  dominated  the  vicinity 
of  the  Coso  geothermal  field,  having  erupted  through  a  Mesozic  basement  highland. 
This  study  focuses  on  the  most  voliuninous  and  areally  extensive  units,  basalts  and 
andesites  of  the  Wild  Horse  Mesa.  These  areally  extensive  units  provide  constraints  on 
fault  offsets  and  fault  kinematics  (vertical-axis  block  rotation). 

2.2.3  Shearing  history 

Wild  Horse  Mesa  provides  a  time-integrated  record  of  shearing  in  the  form  of 
outcrop-scale  brittle  faults.  Wild  Horse  Mesa  is  characterized  by  a  systematic  decrease 
in  elevation  from  east  to  west.  The  mesa  consists  of  relatively  flat,  slightly  west-tilting 
topographic  benches  that  step  down  in  elevation  across  dominantly  west-facing 
normal-fault  scarps.  Canyons  cut  the  mesa  at  high  angles  to  the  overall  topographic 
relief,  forming  relatively  deep,  largely  northeast-trending  scars  in  which  the  basement 
of  the  mesa  is  well  exposed.  Topographically  tiie  lower  reaches  of  these  canyons  cut  the 
western  and  southwestern  parts  of  the  mesa,  whereas  the  upper  reaches  cut  the  eastern 
and  northeastern  parts  of  the  mesa.  Outcrop-scale  faults  have  been  observed  and 
kinematic  constraints  measiured  in  lava  flows  of  (1)  the  lower  section  of  Petroglyph 
Canyon,  (2)  Petroglyph  Canyon  from  near  the  jimction  with  Horse  Canyon  to  its  upper 
reaches  near  flie  end  of  the  "Big  Petroglyphs"  imimproved  road,  (3)  the  lower  section  of 
Dead  End  Canyon,  (4)  the  middle  section  of  Black  Canyon,  and  (5)  the  middle  and 
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upper  sections  of  Sheep  Canyon  (Figure  2.3).  These  data  have  been  exanained  for 
quality  purposes  and  subsequently  divided  on  the  basis  of  location.  Observations  from 
Dead  End  and  Black  Canyons  were  sparse  and  are  not  described  further. 

Data  from  lower  Petrogl)q>h  Canyon,  and  from  Sheep  Canyon  and  the  upper  parts 
of  Petroglyph  Canyon  provide  numerous  kinematic  constraints  on  the  shearing  Idstory 
of  Wild  Horse  Mesa  (Table  2.1  and  Figure  2.4).  These  data  have  been  subdivided  on  the 
basis  of  their  distribution  in  tiie  field.  The  topographically  highest  data  station  in  lower 
Petrogl)q)h  Canyon  is  approximately  300  m  lower  than  the  lowest  data  station  from  the 
upper  Petroglyph  Canyon/Sheep  Canyon  area.  The  data  from  the  region  of  Sheep  and 
upper  Petroglyph  Canyons  have  been  combined  to  form  the  upper  Wild  Horse  Mesa 
data  subset.  The  data  from  fiie  area  aroimd  lower  Petroglyph  Canyon  have  been 
consolidated  into  the  lower  Wild  Horse  Mesa  data  subset. 

The  faults  examined  at  the  upper  and  lower  Wild  Horse  Mesa  display  a  wide  variety 
of  orientations  and  accommodate  normal,  left-lateral,  right-lateral  and  thrust 
displacements  (Figure  2.5).  In  spite  of  the  transtensional  setting,  thrust  faults  make  up  a 
substantial  fraction  of  the  entire  data  set.  It  is  noteworthy  that  the  faults  examined  in 
the  field  are  not  necessarily  spatially  associated  with  the  major  fault  scarps  that  are 
expressed  so  vividly  across  the  step-faulted  terrane.  Field  observations  suggest  that  if 
outcrop-scale  faults  with  readily  interpreted  kinematic  indicators  exist  in  these 
locations,  they  are  largely  covered  with  talus.  High-quality  fault  data  are  most  readily 
available  in  canyon-wall  exposures  immediately  adjacent  to  the  traces  of  major  fault 
scarps  (e.g.,  within  approximately  100  meters  of  where  major  scarps  crossed  canyons). 

2.3  Technical  Approach 

2.3.1  Micropolar  Theory 

Partial  strain  tensors  were  determined  by  treating  the  crust  as  a  micropolar 
continuum  (Twiss  et  al.  1991;  Twiss  et  al.  1993;  Unndi  et  al.  2002).  In  essence  this  entails 
treating  the  crust  as  a  granular  material  in  which  grains  are  free  to  rotate  at  different 
rates  than  the  large-scale  continuum.  The  resulting  model  consists  of  five  parameters, 
three  Euler  angles  that  define  the  orientations  of  the  principal  strain  rates  (i.e.,  > 

di  ^  d^;  lengthening  is  reckoned  positive);  a  deformation-rate  parameter  D  s  (dj  -  d^)/ (d^ 
-  dj),  which  defines  the  relative  magnitudes  of  the  principal  strain  rates;  and  a  relative 
vorticity  parameter  W  =  (CO13  -  Wi3)/0.5(di  -  d^,  which  defines  a  normalized  difference 
between  the  local  vorticity  of  individual  fault  blocks  (O13  and  file  large-scale  vorticity  of 
the  deforming  continuum  W13  about  the  d^  axis. 

This  approach  can  be  used  to  examine  earthquake  focal  mechanisms  (as  in  Figure 
2.1,  and  in  Lewis  et  al.  2003;  Unruh  et  al.  2003;  Unruh  et  al.  2002;  Unruh  et  al.  1996)  or 
fault-slip  data  (as  in  Figures  3.3-3.6,  and  in  Unruh  &  Twiss  1998).  The  fundamental 
difference  in  these  two  applications  is  (1)  nodal  plane  ambiguity  in  the  former  and  no 
such  ambiguity  in  the  latter,  and  (2)  modeling  instantaneous  deformation  in  the  former 
(i.e.,  essentially  strain  rate)  and  modeling  part  of  the  time-integrated  (i.e.,  finite)  strain 
in  the  latter. 

The  robustness  of  the  best-fitting  strain  tensors  is  evaluated  using  bootstrap  statistics 
to  constrain  the  variation  among  the  five  model  parameters.  D  values  statistically  less 
than  or  greater  than  0.5  indicate  deviatoric  transtension  or  deviatoric  transpression, 
respectively.  W  values  statistically  different  than  zero  indicate  crustal  block  rotation 
fiiat  is  different  than  that  of  the  large-scale  continuum.  For  seismogenic  deformation  W 
is  interpreted  in  terms  of  rotation  rates  because  the  seismic  record  essentially  provides 


Purchase  Request  Number:  N68936-01-C-0094 


7 


constraints  on  incremental  deformation.  For  the  rock  record,  W  is  considered  to  reflect 
differences  in  finite  block  rotation  relative  to  that  of  the  large-scale  continuum. 

The  seismogenic  deformation  at  Wild  Horse  Mesa  was  evaluated  starting  with  the 
results  of  Unruh  et  al.  (2002).  Two  primary  analyses  were  performed:  1)  tests  of 
differing  deptiis  for  the  transition  from  crustal  thinning  at  shallow  levels  to  horizontal 
shearing  at  deeper  levels,  and  2)  statistical  evaluations  of  non-plane  strain  and  non-zero 
relative  vorticity. 

In  determining  best-fitting  strain  tensors  to  fault-slip  data  the  data  were  sorted  in  a 
variety  of  ways  in  order  to  maximize  the  possibility  that  the  data  reflected 
homogeneous  strain.  Given  the  time-integrated  nature  of  the  fault  data,  this  sorting 
requires  multiple  approaches.  Four  primary  means  of  subdividing  the  data  have  been 
used.  First  the  data  are  examined  on  the  basis  of  where  they  occur  in  the  field,  that  is, 
spatially.  Due  to  the  sparseness  of  kinematic  indicators  this  approach  is  somewhat 
limited.  Nonetheless,  two  stations  (203  and  213,  see  Table  2.1)  yielded  enough  fault-slip 
data  to  warrant  this  approach.  Second,  pooled  data  are  examined  in  hopes  of  finding 
multiple  strain  geometries  that  are  viable.  These  geometries  are  then  used  to  subdivide 
that  data  based  on  the  affinity  of  a  data  point  for  one  strain  tensor  over  another  using  a 
minimum  misfit  criterion,  lids  approach  is  "blind  to"  the  spatial  distribution  of  the 
fault  data.  Third,  fault  data  are  subdivided  on  the  basis  of  tlie  style  of  deformation 
accommodated.  In  this  approach,  for  example,  normal  faults  and  strike-slip  faults  are 
separated.  Fourth,  clusters  of  similarly  oriented  kinematic  P  and  T  axes  are  used  to 
subdivide  the  data.  This  is  typically  done  using  Kamb  contour  plots  of  the  P  and  T 
axes. 

The  primary  spatial  subdivision  of  the  brittle  fault  data  is  based  on  position  within 
Wild  Horse  Mesa.  The  topographically  highest  observation  from  lower  Wild  Horse 
Mesa  is  approximately  300  m  lower  than  d\e  lowest  observation  from  upper  Wild  Horse 
Mesa.  Moreover,  the  fault  traces  mapped  by  Duffield  et  al.  (1980)  display  distinct 
^  patterns  in  these  two  areas.  Lower  Wild  Horse  Mesa  as  we  have  defined  it  is 

dominated  by  north  to  slightly  northeast  trending  fault  traces.  Upper  Wild  Horse  Mesa 
is  likewise  characterized  by  this  fault  trace  geometry,  but  also  includes  a  set  of 
nothwest-trending  fault  traces  (Figiure  2.3).  The  fault  data  from  lower  Wild  Horse  Mesa 
are  heterogeneous  and  only  one  homogenous  subset  is  apparent,  as  discussed  below. 

The  data  from  upper  Wild  Horse  Mesa  have  been  divided  using  a  combination  of 
approaches.  The  fault  slip  data  from  station  203  suggest  two  primary  styles  of 
deformation,  crustal  thickening  and  horizontal  shearing.  These  two  styles  of  strain  are 
also  apparent  in  the  data  from  other  parts  of  upper  Wild  Horse  Mesa.  Because  stations 
203  and  213  are  within  about  1  km  of  one  another  and  within  about  2  km  of  all  of  the 
remaining  observations  from  upper  Wild  Horse  Mesa,  it  is  considered  reasonable  to 
initially  pool  all  of  the  data  from  the  upper  Wild  Horse  Mesa.  These  data  were  then 
subdivided  on  the  basis  of  the  type  of  faulting,  using  dip-slip  dominated  normal 
faulting,  strike-slip  dominated  horizontal  shearing  and  dip-slip  dominated  thrust 
faulting  as  catagories. 

2.3.2  X-Ray  Fluorescence 

The  intention  of  this  task  is  to  fingerprint  eruptive  units  in  order  to  compare 
paleomagnetic  remanence  directions  from  the  same  unit  at  different  localities. 
Separation  of  imits  by  appearance  in  the  field  or  in  thin  section  is  not  sufficiently 
diagnostic  or  rigorous.  Therefore,  we  sought  geochemical  means. 
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Throughout  this  report  the  term  eruptive  unit  signifies  a  lava  flow  or  flows 
exhibiting  a  (within-locality)  uniform  paleomagnetic  remanence  direction  and,  usually, 
a  tightly  defined  chemistry.  These  criteria  mdicate  rapid  eruption  and  allow  rigorous 
correlation.  Identification  of  units  permits  the  measurement  of  relative  vertical-axis 
rotation  of  fault  bounded  blocks  between  localities.  In  the  process  a  much  more 
detailed  imderstanding  of  stratigraphic  relations  has  emerged  than  previous  studies 
have  yielded  (Novak  and  Bacon,  1986;  Duffield  and  Bacon,  1981).  Our  data  indicate 
that  previous  models  of  the  geochemical  evolution  of  the  Pliocene  volcanic  field  at  Wild 
Horse  Mesa  require  amendment. 

Using  the  method  of  Johnson  et  al.  (1999)  we  prepared  65  lava  samples  for  major 
(oxides  of  Na,  Mg,  Al,  Si,  P,  K,  Ca,  Ti,  Mn,  and  Fe)  and  trace  (Rb,  Sr,  Y,  Zr,  and  Nb) 
element  analysis  by  X-Ray  Fluorescence  (XRF)  spectroscopy.  Table  2.2  shows  raw 
results  for  these  samples.  Supplementary  data  on  V,  Ni  and  Cr  for  selected  samples  are 
included  in  Appendix  A.  Each  sample  was  processed  and  analyzed  in  duplicate  to 
assess  reproducibility  of  preparation  procedures.  The  standard  deviation  reported  is 
that  calculated  from  measurement  of  the  duplicate  samples  plus  any  repeated 
measurements.  The  standard  deviation  of  individual  measurements  should  not  be 
confused  with  the  true  chemical  variability  throughout  a  flow  or  eruptive  imit. 

Major  elements  and  trace  elements  were  analyzed  using  2  different  settings  on  the 
same  Rh  x-ray  tube  (45  kV,  50  mA  for  major  elements  and  90  kV,  30  mA  for  trace 
elements).  During  analysis  of  a  sample  batch,  several  quality  assurance/ quality  control 
samples  were  nm.  International  standards  BIR-1  and  AGV-1  were  each  analyzed  as  a 
calibration  check  once  or  twice  during  each  3-day  analysis  session.  Counting  standard 
RLGM  was  analyzed  at  least  once  each  day.  Remeasurement  of  selected  samples  from 
previous  batches  further  ensured  consistency  between  sessions.  These  steps  were  all 
taken  to  ensure  die  accuracy  of  XRF  analysis. 

To  gauge  the  real  chemical  variation  of  flows  within  eruptive  units,  we  typically 
selected  two  samples  from  each  eruptive  unit  at  a  given  locality  for  analysis.  Therefore, 
the  spread  on  concentrations  for  a  given  unit  reflects  the  real  chemical  variation,  both 
due  to  initial  emplacement  conditions  and  variable  (though  usually  minimal) 
hydrothermal  alteration  and  weathering.  Element  concentrations  were  normalized  to 
100%,  ignoring  volatile  content  Johnson  et  al.,  1999),  usually  an  indicator  of  post¬ 
emplacement  hydrothermal  and  weathering  processes  rather  than  parent  magma 
content. 

Geochemistry  alone  is  insufficient  to  fully  distinguish  eruptive  units.  Some  lava 
flows  with  opposite  paleomagnetic  polarity  exhibit  similar  geochemistry.  Thus,  we 
fingerprint  units  using  combined  geochemistry  and  paleomagnetism.  Details  of 
paleomagnetic  techniques  and  results  wiU  be  discussed  in  subsequent  sections. 

Some  workers  have  used  geochemistry  of  major  elements  (e.g.,  Monastero  et  al., 
1997)  sometimes  combined  with  paleomagnetism  (e.g..  Wells  et  al.,  1989)  to  correlate 
widely  separated  outcrops  of  the  same  lavas.  However,  we  prefer  correlation  using 
mainly  trace  elements  (Rb,  Sr,  Y,  Zr,  and  Nb).  This  is  because  major  element 
geochemistry  can  vary  significantly  during  the  coiurse  of  an  eruption  due  to 
differentiation  processes  and  zoned  magma  chambers.  Since  we  require  separation  of 
units  with  similar  bulk  chemistry  (and  appearance)  tibieir  bulk  chemistries  overlap.  For 
this  reason,  we  move  to  trace  elements.  Partition  coefficients  into  basaltic  minerals  are 
low  in  these  elements  and  of  similar  magnitude,  such  that  their  concentration  and 
especially  their  ratio  during  an  individual  eruption  is  likely  to  be  invariant  for  basalts 
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(and  increase  in  variability  with  larger  silica  content)  even  in  the  face  of  some 
^ferentation. 


2.3.3  Paleomagnetism 

One  method  for  revealing  the  kinematics  of  fault  blocks  in  the  study  area  is  to 
measure  finite  rotations  using  paleomagnetic  techniques.  Two  methods  of  measuring 
block  rotation  are  described.  Raw  paleomagnetic  data  are  not  included  here  due  to 


space  constraints.  These  data  can  be  foimd  in  Pluhar's  thesis  (Pluhar,  2003) 
Paleomagnetism  of  Sediments 

In  the  most  simple  case,  paleomagnetic  measurements  of  vertical-axis  rotations  are 
made  on  sediments.  Sampling  through  a  stratigraphic  section  typically  allows 
averaging  of  secular  variation  of  Earl's  magnetic  field.  In  addition  there  is  no 
ambiguity  on  the  exact  bedding  correction  to  apply,  since  all  but  fanglomerate  is 
deposited  close  to  horizontal.  Thereby,  the  mean  direction  of  a  group  of  samples 
through  such  a  section  points  to  paleogeographic  north  for  the  tectonic  block  on  which 


the  sediments  ride. 


We  sampled  sediments  at  three  localities  (see  Figure  2.2  for  locations).  The  Haiwee 
Reservoir  (HR)  lacustrine  section  yielded  poor  directional  results  and  will  be  discussed 
no  further.  The  White  Hills  section  B  (WHB)  is  also  lacustrine  and  thought  to  be 
Pleistocene  (Duffield  et  al.,  1980).  BoA  of  these  lacustrine  sections  were  sampled  by 
collection  of  large,  oriented  handsamples  which  were  then  subdivided  into  multiple 
specimens.  The  Coso  Wash  (CWA)  section  spans  part  of  the  Coso  formation 
fanglomerate  member  and  is  thought  to  be  late  Pliocene  in  age  (Duffield  et  al.,  1980). 
Each  sample  collected  from  this  locality  was  individually  carved,  a  plastic  box  placed 
over  the  in  situ  sample,  and  orientation  recorded.  This  locality  is  unique  in  that  it 
contains  sediment  as  the  majority  of  section  beneath  an  angular  unconformity  capped 
by  pyroclastics  and  a  dacite  (Duffield  and  Bacon,  1981,  mapping  unit  Tdo).  This  allows 
a  direct  radiometric  tie  point,  discussed  in  the  geochronology  section,  to  the  magnetic 
reversal  timescale.  For  the  capping  dacite,  paleomagnetic  techniques  described  in  the 
next  section  for  lavas  were  applied. 

Oriented  sediment  samples  were  analyzed  by  standard  paleomagnetic  techniques 
including  tiiermal  and/or  alternating  field  (AF)  stepwise  demagnetization  experiments 
and  principal  component  analysis  (PCA)  (IQrschvi^,  1980;  Cogne,  2003)  of  the 
resulting  demagnetization  padis.  For  sediments,  thermal  demagnetization  was  the 
preferred  technique  because  of  the  possibility  of  secondary  chemical  remanent 
magnetization  components.  PCA  yields  the  primary  or  characteristic  remanent 
magnetization  (ChRM)  direction  for  each  specimen.  Where  multiple  specimens  were 
derived  from  a  single  oriented  hand  sample,  all  ChRMs  determined  for  the  specimens 
were  averaged  to  calculate  a  mean  ChRM  for  the  sample.  For  each  locality,  ihe  mean  of 
an  independently-oriented  ChRMs  yields  paleogeographic  north. 

For  locality  WHB,  standard  tilt  corrections  were  applied  to  restore  beds  to  original 
horizontality.  However,  for  CWA,  the  fanglomerate  locality,  the  assumption  of  original 
horizontahty  is  erroneous.  Alluvial  fans  typically  exhibit  slopes  (i.e.,  original  dips)  of 
5°-10°  as  do  the  modem  fans  in  Coso  Wash.  Thus  at  CWA  we  restore  beds  tilting  25° 
southeastward  to  either  5°  or  10°  southeastward,  the  minimum  restoration  required  that 
is  still  geologically  reasonable  given  no  change  in  paleo-drainage  direction.  Restoration 
to  original  dips  in  some  other  direction  for  the  most  part  result  in  even  larger  clockwise 
rotations  than  those  observed  for  the  minimum  corrections  (see  results  section). 
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Paleomagnetism  of  Lavas 

The  major  Cenozoic  rock  types  available  in  the  study  area  are  basalt  and  andesite. 
Lavas  are  often  less  than  ideal  for  averaging  secular  geomagnetic  variation  because 
many  flows  can  be  emplaced  in  a  very  short  period,  yielding  a  thick  stratigraphic 
section  representing  the  same  geologic  instant,  and  tiierefore,  paleomagnetic  direction. 
However,  we  utiUze  such  rocl«  because  they  present  other  advantages  for 
paleomagnetic  studies.  Lavas  very  faithfully  record  Earth's  magnetic  field,  and 
generally  this  paleomagnetic  signal  can  be  measured  more  accurately  and  precisely 
than  tiiat  in  sediments.  Therefore,  if  a  laterally-extensive  lava  flow  can  be  identified, 
measuring  the  directional  differences  from  one  locality  to  another  reveals  the  relative 
rotation  between  localities  to  wifiiin  about  ±  5®. 

In  this  discussion  we  refer  to  individual  cooling  units  as  lava  flows,  and  groups  of 
lava  flows  which  exhibit  strong  evidence  of  rapid  emplacement  as  eruptive  units  (called 
flow  fields  by  some  workers).  Evidence  of  rapid  emplacement  are:  statistically 
indistinguishable  remanence  direction  between  successive  flows  and/or  identical 
geochemical  "fingerprint"  (see  preceding  section). 

We  sampled  22  localities  of  late  Pliocene  and  early  Pleistocene  lavas.  Localities 
varied  in  character,  representing  anywhere  from  one  to  twenty-one  individual  lava 
flows.  At  most  localities  (BC,  CWA,  UDC,  LB,  PC,  HS,  HSB,  HSD,  SP,  CF,  CHS,  RCB, 
RCC,  RCD,  RCE)  six  or  more  oriented  drill  cores  were  collected  from  each  lava  flow.  At 
a  few  localities  (AL,  RC,  BD)  as  few  as  three  samples  per  flow  were  collected  due  to 
time  constraints  or  logistical  limitations.  At  one  locality  we  collected  oriented  hand 
samples  for  recormaissance  purposes  (LDC).  These  data  are  identified  as  such  in 
subsequent  sections.  Results  from  localities  CMA,  CMC,  and  MSA,  yielded  data  of 
minor  importance  which  will  be  discussed  in  passing  only.  For  all  lava  samples  except 
handsamples,  sun  compass  or  sight  point  corrections  were  applied  to  accoxmt  for  local 
magnetic  anomalies. 

Lava  specimens  drawn  from  the  samples  were  stepwise  tibermal  or  AF 
demagnetized.  Specimens  typically  exhibited  imivectorial  decay  to  the  origin  imless 
affected  by  lightning  strikes.  Since  AF  demagnetization  is  far  more  effective  than 
thermal  for  isolating  magnetization  components  in  lightning-affected  samples,  it  was 
the  favored  method  for  lavas.  Demagnetization  data  were  analyzed  using  PCA,  and 
mean  directions  were  calculated  for  each  lava  flow  from  ti\e  sample  ChRMs  for  that 
flow.  For  eruptive  imits  that  contain  multiple  lava  flows  and  for  which  flow  means  are 
statistically  indistinguishable  from  one  another,  the  flow  mean  directions  were 
averaged  to  yield  a  unit  mean  direction  for  that  locality.  Thus  for  the  tectonic  analysis 
we  compare  unit  mean  directions  between  localities. 

Attitude  of  lava  flow  tops  were  recorded  as  an  estimate  for  paleohorizontal  at  flow 
emplacement.  These  attitudes  were  measured  in  the  field  where  possible  or  taken  from 
aerial  photographs  using  a  parallax  bar  by  standard  mettiods.  However,  the 
assumption  of  original  horizontality  of  flow  tops  is  obviously  a  simplification,  since 
lava  flows  downslope.  For  tilt  corrections  in  this  study  we  restore  flow  tops  to 
horizontal  unless  otiierwise  noted. 

In  addition  to  lavas,  we  collected  and  anlayzed  airfall  piunice  samples  from  two 
localities  (RHT  and  RDT).  These  samples  appear  to  have  acquired  much  of  their 
magnetization  during  flight  or  before  achieving  a  final  position  in  the  agrading  pumice 
pile  and  thus  carry  many  components  and  appear  imstable  at  moderate  to  high 
demagnetization  steps.  In  addition,  thermal  demagnetization  tends  to  break  samples 
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apart  by  differential  expansion  of  the  large  micas  contained  therein.  Thus,  these  low 
quality  data  will  be  discussed  no  further. 

2.3.4  Geochronology 

Radiometric  dating  was  undertaken  in  order  to  better  constrain  the  time  spanned  by 
the  finite  rotation  determined  paleomagenticaUy  and  by  fault-slip  inversion.  Sebastien 
Nomade  imder  the  supervision  of  Paul  Renne  at  Berkeley  Geochronology  Center  (BGG) 
conducted  6  analyses  for  ^Ar/^’Ar  geochronology  in  accordance  with  the  statement  of 
work.  Techniques  and  many  of  the  details  of  this  effort  can  be  found  in  BGC's  attached 
report  (Appendix  B).  Also  presented  are  data  for  eight  additional  analyses  conducted 
by  Jonathan  M.G.  Glen  and  Sebastien  Nomade  at  BGC  (raw  data.  Appendix  B)  which 
have  not  yet  been  published  elsewhere. 

In  conducting  the  six  analyses  in  satisfaction  of  our  contract,  it  was  foimd  that 
quartz  and  perhaps  plagioclase  xenocryst  contamination  of  these  samples  was 
widespread  in  the  Pliocene  lavas.  For  this  reason  groundmass  separates  were  favored. 
This  fact  should  be  kept  in  mind  when  judging  the  reliability  of  previous  K/ Ar  whole 
rock  ages  (e.g.,  Duffield  et  al.,  1980). 

3.0  Results 

3.1  Seismogenic  Strain 

The  seismogenic  deformation  at  Wild  Horse  Mesa  was  evaluated  starting  with  the 
results  of  Unrvdi  et  al.  (2002).  Two  primary  analyses  were  performed:  (1)  tests  of 
differing  depths  for  the  transition  from  crustal  thinning  at  shallow  levels  to  horizontal 
shearing  at  deeper  levels,  and  (2)  statistical  evaluations  of  non-plane  strain  and  non¬ 
zero  relative  vorticity.  The  results  of  the  first  analyses  indicate  that  5  km  indeed 
represents  the  transition  from  crustal  thinning  to  horizontal  shearing.  The  results  of  the 
second  analyses  reveal  that  seismogenic  deformation  reflects  plane  strain  partitioned  in 
depth  (i.e.,  D  =  0.5)  and  that  block  rotation  occurs  at  rates  different  than  those  for  the 
large-scale  continuum  (i.e.,  W  is  non-zero).  For  the  shallow  events  (Figure  3.1),  W  = 
0.456  (+0.283/-0.232,  at  95%  confidence  here  and  throughout,  except  if  noted  otherwise) 
and  for  the  deeper  events  (Figure  3.2),  W  =  0.647  (+0.200/-0.426).  For  deformation  at 
shallow  levels  tins  value  is  consistent  with  north-viewed  coimterclockwise  rotation,  or 
clockwise  rotation  at  a  rate  slower  than  the  macrovorticity.  Given  the  geometry  of  the 
well-developed  west-facing  fault  scarps  across  Wild  Horse  Mesa,  clockwise  block 
rotation  is  the  preferred  interpretation.  For  deformation  at  deeper  levels,  the  W  value  is 
consistent  with  coimterclockwise  vertical-axis  rotation,  or  clockwise  vertical-axis 
rotation  at  a  rate  slower  than  the  macrovorticity.  Considering  that  the  seismic  events 
define  a  northwest  trending  planar  source  zone,  and  that  the  boundary  conditions  are 
dominated  by  dextral  shear,  the  latter  interpretation  is  preferred. 

Seismogenic  strain  at  both  structural  levels  is  consistent  with  expectations  given  the 
transtentional  setting.  For  shallower  and  deeper  seismdty,  dj  is  subhorizontal  and  west 
soutiiwest-east  northeast  trending.  Crustal  thinning  at  shallow  levels  is  characterized 
by  a  subvertical  dj,  whereas  horizontal  dextral  shearing  at  deep  levels  is  characterized 
by  subvertical  dj. 

3.2  Exposed  Outcrop-Scale  Brittle  Faults 

The  lower  Wild  Horse  Mesa  faults  that  we  examined  are  characterized  by 
heterogeneously  distributed  kinematic  P  and  T  axes  as  well  as  diverse  fault  plane  and 
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lineation  orientations  (Figure  2.5):  To  date,  our  attempts  to  parse  these  data  into 
homogeneous  subsets  have  yielded  limited  success.  Inversion  of  the  collective  fault 
data  for  which  we  have  high  confidence  yields  an  oblique  strain  tensor  that  indicates 
shallow  south-plimging  extension  and  shallow  west-plunging  shortening  (not  shown). 
The  south-trending  maximxim  stretching  direction  (d^)  is  quite  different  tiian  the 
subhorizontal,  nearly  east-west  trend  that  is  characteristic  of  dj  for  seismogenic 
deformation  in  the  region  (Figures  3.1  and  3.2).  Nonetheless,  the  dominant  solution 
geometries  suggest  that  the  shallow  west  northwest  pltmging  P  axis  cluster  apparent  in 
Figure  2.5  reflects  the  existence  of  a  homogeneous  subset  of  fault  slip  data.  Using  this  P 
axis  maximum  as  a  means  of  subdividing  the  data  set  results  in  the  isolation  of  a  fairly 
well-behaved  data  subset  (Figure  3.3).  In  spite  of  the  fact  that  this  subset  appears 
statistically  homogeneous  (i.e.,  the  final  positions  for  principal  axes  are  clustered  in 
Figure  3.3c),  the  resulting  strain  tensor  is  difficult  to  interpret  in  tire  context  of  the 
regional  tectonics. 

The  upper  Wild  Horse  Mesa  faults  display  a  less  heterogeneous  distribution  of 
kinematic  P  and  T  axes  than  do  the  faults  of  lower  Wild  Horse  Mesa.  Moreover,  the 
fault  planes  are  somewhat  organized  with  two  sets  apparent.  The  first  set  is  dominated 
by  subvertical  faults  that  are  either  west  or  northwest  striking.  The  second  set  is 
dominated  by  faults  that  stike  north  and  dip  moderately  eastward.  Nonetheless,  these 
faults  accommodate  crustal  thinning,  crustal  thickening  and  horizontal  shearing. 
Generally  speaking,  these  faults  display  lineations  tiiat  are  either  approximately  down 
dip  or  approximately  strike  parallel.  As  a  consequence  the  faults  have  been  subdivided 
into  subsets  on  the  basis  of  style;  that  is,  normal  faults,  thrust  faults  and  strike-slip 
faults. 

Inverse  modeling  of  the  faults  for  which  we  have  high  confidence  in  our  field 
observations  yields  three  strain  tensors.  The  results  for  the  normal  faults  (Figure  3.4) 
and  for  the  thrust  faults  (Figure  3.5)  are  statisically  less  robust  than  for  the  strike-slip 
faults  (Figure  3.6)  owing  to  die  smaller  niunbers  of  observations.  As  a  consequence,  it  is 
not  possible  to  constrain  non-plane  strain  or  non-zero  relative  vorticity.  In  contrast,  at 
95%  confidence  the  strike-slip  faults  suggest  a  strong  bias  toward  a  prolate  strain  tensor 
and  non-zero  relative  vorticity  (Figure  3.6,  c  and  e).  Accepting  90%  confidence,  it  can  be 
argued  that  these  faults  indeed  accommodate  prolate  strain,  consistent  with 
trai\stension,  and  non-zero  relative  vorticity,  consistent  with  coimterclockwise  vertical- 
axis  rotation,  or  clockwise  vertical-axis  rotation  at  a  rate  slower  than  the  macrovorticity. 
Given  the  regional  history  of  dextral  transtension,  the  latter  interpretation  of  the  relative 
vorticity  result  is  preferred. 

The  inversion  results  for  the  upper  Wild  Horse  Mesa  faults  are  internally  consistent 
and  in  general  agreement  with  residts  for  contemporary  strain.  The  orientation  of  dj  is 
consistent  for  the  strike-slip  and  normal  faults;  namely  subhorizontal  and  west 
southwest-east  northeast  trending  (compare  Fi^es  3.4  and  3.6).  Similarly,  for  the 
strike-slip  faults  and  the  thrust  faults,  dj  is  consistently  subhorizontal  and  north-south 
trending  (compare  Figures  3.5  and  3.6).  Although  these  three  fault  subsets  indicate 
strain  that  has  partitioned  at  sub-km  scales,  the  overall  geometries  of  strain  are  in 
accord  with  the  seismogenic  strain  (compare  Figures  3.1  and  3.2  with  Figures  3.4-3.6). 
The  oidy  notewortihy  exception  is  the  thrust  faults,  which  reveal  crustal  thickening  that 
is  not  apparent  in  the  seismogenic  deformation.  The  existence  of  crustal  thickening, 
however,  is  not  unexpected  because  of  kinematic  compatibility  requirements  during 
vertical-axis  block  rotation  (Dewey  2002). 
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3.3  Unit  Correlation  from  X-Ray  Fluorescence  Geochemist^ 

In  order  to  map  the  pattern  of  block  rotation  it  is  necessary  first  to  define  laterally 
extensive  eruptive  units  on  which  to  perform  paleomagnetic  analyses.  Duffield  and 
Bacon  (1981)  defined  much  of  the  greater  Wild  Horse  Mesa  as  a  single  unit  called  Basalt 
of  Petroglyph  Canyon.  We  subdivide  this  gross  compositional  unit  into  multiple  new 
eruptive  units  each  spanning  a  period  shorter  than  that  necessary  for  noticeable  secular 
variation  of  Earth's  magnetic  field  (about  200  years)  and  sharing  a  similar  geochemistry. 
Each  newly  defined  uiut  probably  represents  the  outpouring  from  a  single  vent  during 
a  geologically  short-lived  eruption. 

As  stated  previously,  trace  element  concentrations  are  more  likely  to  be  invariant 
during  an  eruption  than  major  elements,  because  they  are  little  affected  by 
differentiation.  Figure  3.7  depicts  Sr/Zr  versus  Rb/Nb  ratios  combined  with 
paleomagnetic  polarity  for  Pliocene  Wild  Horse  Mesa  volcanics.  These  ratios  are 
chosen  because,  of  all  combinations  of  the  five  trace  elements  analyzed,  they  provide 
the  best  separation  of  units  into  discrete  fields.  This  separation  into  distinct  fields 
suggests  a  unique  chemical  fingerprint  for  each.  Magnetic  polarity  is  indicated  on  the 
figures  as  open  symbols  for  reversed  polarity,  solid  symbols  for  normal,  and  characters 
(asterisk,  cross,  etc.)  for  unknown  polarity  (e.g.  for  samples  composed  of  drillhole 
cuttings).  Using  the  diagram,  we  assign  new  eruptive  unit  names  based  on  the 
probable  eruptive  soirrce  area  of  each  unit.  Classification  into  basalt,  andesite,  and 
dacite  anticipates  major  element  data  presented  in  Figures  3.9-3.22. 

We  subdivide  the  former  Basalt  of  Petroglyph  Canyon  into  a  redefined  Petroglyph 
Canyon  basalt.  Silver  Peak  basalt.  Cole  Flat  basalt,  Haiwee  Spring  basalt,  and  Black 
Canyon  flows  (some  of  which  were  liunped  into  Duffield  and  Bacon's  Petroglyph 
Canyon  imit).  Petroglyph  Canyon  basalt  is  easily  distinguished  from  Black  Canyon 
flows  by  polarity  but  not  geochemistry.  On  the  other  hand,  Petroglyph  Canyon  basalt 
is  separated  from  Silver  Peak  basalt  by  trace  element  chemistry  but  not 
paleomagnetism.  As  this  example  shows,  both  geochemistry  and  paleomagnetism  are 
required  to  separate  the  units  studied  here. 

It  should  be  noted  that  within  the  Black  Canyon  unit,  a  dacite  intercalated  (and  not 
intruded)  amongst  11  basalts,  while  visually  and  geochemically  distinguishable,  shares 
the  same  imusual  remanence  direction  (see  section  3.4).  Accordingly  it  could  be 
separated  as  a  distinct  stratigraphic  imit.  Due  to  its  localized  nature  and 
contemporaneity  with  flows  above  and  below  based  on  remanence,  we  refrain  from 
doing  so. 

Another  important  result  is  the  separation  of  the  Coso  Basin  andesite  (previously 
dubbed  Basalt  of  Coso  Basin  by  Whitmarsh,  1996)  from  any  other  imit  sampled,  using 
both  geochemistry  and  paleomagnetic  remanence  direction.  Field  evidence  also 
suggests  that  ihe  only  surficial  locality  of  this  eruptive  unit  is  also  the  source  region 
(local  dikes  and  intrusions  of  the  same  polarity  are  foimd  at  the  locality). 

In  addition,  geochemistry  suggests  correlation  of  lavas  penetrated  at  1300'  and  1000' 
below  grade  by  drill  hole  35-2  (Figure  3.7)  with  Black  Canyon  flows  and  Petroglyph 
Canyon  basalt  respectively.  This  tentative  correlation  is  discussed  further,  below. 

In  the  Renegade  Canyon  area  we  analyzed  lavas  designated  Basalt  of  Renegade 
Canyon  (QTbr)  and  divided  into  subunits  by  Duffield  and  Bacon  (1981).  Figure  3.8 
shows  trace  element  data  for  lavas  of  this  area,  and  refines  ihe  previous  subimit 
definition  geochemically,  allowing  widespread  correlation  of  these  units  across  the  area. 
Most  of  these  units  are  reversed  polarity  with  very  similar  remanence  directions  and 
tihus  they  are  defined  largely  on  the  basis  of  geochemistry.  Results  from  Renegade 
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Canyon  suggest  that  the  geochemistry  of  eruptive  units  tiiere  varies  significantly  within 
eruptive  units  and  even  wilhin  individual  lava  flows.  For  example,  samples  of  QTbrg 
from  a  single  locality  and  QTbr4  from  two  different  localities  vary  significantly.  This 
introduces  some  uncertainty  into  correlations  around  the  Renegade  Canyon.  Question 
marks  enclosing  QTbri  on  Figure  3.8  indicate  the  imcertainty  due  to  trace  and  major 
element  geochemical  variation  as  well  as  paleomagnetic  results. 

Major  element  and  trace  element  geochemistry  are  plotted  versus  MgO  on  Figures 
3.9-3.22  and  strongly  support  units  defined  for  Wild  Horse  Mesa  lavas  (Figure  3.7)  by 
tile  tight  grouping  of  data  points  within  each  eruptive  unit.  These  data  also  support 
units  QTbrz  and  QTbr4  as  defined  (Figure  3.8).  However,  QTbri  is  a  questionable 
grouping  because  of  the  lack  of  paleomagnetic  data  and  the  spread  on  major  element 
geochemistry  for  locality  RC  lava  flow  3  and  the  bottom  flows,  compared  to  the  QTbri 
type  locality.  As  will  be  shown  in  section  3.4  this  uncertainty  is  immaterial  because 
there  appears  to  be  no  measurable  relative  rotation  among  localities  of  Renegade 
Canyon  over  the  age  of  tiie  lavas. 

Results  from  major  element  analysis  also  call  into  question  possible  correlations  of 
cuttings  samples  from  drill  hole  35-2  from  which  no  paleomagnetic  analyses  are 
possible.  These  geochemical  data,  for  the  most  part,  fall  far  off  the  trend  of  data  from  all 
surface  samples.  The  proximity  of  this  drill  hole  to  the  geothermal  area  and  the 
detection  in  nearby  drillholes  of  geothermal  waters  suggests  that  the  xmusual 
geochemistry  of  these  lavas  results  from  hydrothermal  alteration.  Such  alteration  is 
known  to  strongly  deplete  samples  of  Mg  at  mid-ocean  ridges.  Therefore,  while 
geothermal  water  appears  to  have  depleted  these  units  in  Mg  and  probably  other 
mobile  major  elements,  we  suggest  that  their  trace  element  concentrations  may  be 
largely  imperturbed.  Furthermore,  the  ratio  of  trace  elements  analyzed  would  be  even 
more  stable  because  of  their  similar  chemical  behavior.  Ergo,  we  retain  correlation  of 
tile  1000'  and  1300'  samples  of  hole  35-2  with  Petroglyph  Canyon  basalt  and  Black 
Canyon  flows. 

The  lateral  extent  of  major  new  or  redefined  units  described  here  are  depicted  in 
map  view  in  Appendix  C. 

3.4  Paleomagnetism 

3.4.1  Sedimentary  Sections 

Sedimentary  localities  offer  the  opportunity  to  measure  vertical  axis  rotation  relative 
to  paleo-geographic  north,  i.e.  absolute  rotation.  Figure  3.23  shows  data  for  localities 
ON  A  and  WHB.  The  dispersion  of  individual  sample  ChRMs  about  the  mean  is  the 
result  of  secular  variation  of  Earth's  magnetic  field  plus  any  infidelity  in  remanence 
recording  by  the  sediments. 

Reliable  data  from  the  WHB  (Figure  3.23A)  are  all  of  reversed  polarity,  consistent 
with  previous  age  constraints  placing  it  in  the  early  Pleistocene  (Duffield  et  al.,  1980). 
The  tilt-corrected  mean  direction  of  D  (declination)  =  175.5°,  I  (inclination)  =  -21.4°  and 
Otgs  (region  that  is  95%  likely  to  enclose  the  true  mean  direction)  =10.1°  for  13 
independently  oriented  handsamples  (N=13)  is  statistically  indistinguishable  from 
geographic  south  but  is  significantly  shallowed  in  inclination.  This  shallowing  from  an 
expected  mean  inclination  of  -55.5°  is  likely  due  to  two  factors;  1)  shallowing  due  to  ti\e 
tendancy  of  magnetized  tabular  hematite,  which  is  ordinarily  magnetized  in  the  plane 
of  the  grain,  to  lie  flat.  Hematite  is  the  main  remanence  carrier  in  the  most  shallowed 
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samples.  2)  post-sedimentation  compaction  as  corroborated  by  a  horizontally-foliated 
fabric  in  anisotropy  of  magnetic  susceptibility  measurements  of  these  sediments. 

The  CWA  section  is  entirely  normal  polarity,  consistent  with  the  inferred  age  of 
Coso  Wash  fanglomerate  of  around  3  Ma  (Duffield  et  al,  1980).  Figure  3.23b  depicts 
ChRMs  for  independently-oriented  fanglomerate  samples,  a  mean  ChRM  for  multiple 
specimens  from  an  oriented  handsample  of  pyroclastics  immediately  beneath  the  dadte 
and  the  mean  direction  of  the  capping  dadte  based  on  four  independently-oriented 
(and  averaged  prior  to  inclusion  into  a  locality  mean)  drill  core.  We  apply  two  possible 
tilt  corrections  in  Figure  3.23B,  restoring  the  fanglomerate  to  original  dips  of  10°  or  5°. 
We  err  on  the  conservative  side  (less  correction)  and  choose  the  former  as  the  best  tilt 
correction.  The  proximity  of  Mesozoic  basement  suggests  the  Pliocene  moimtain  front 
producing  this  fanglomerate  was  similarly  proximal  and  thus  the  steep  original  dip  is 
warranted. 

Thus,  the  CWA  locality  mean  direction  is  D  =  21.6°,  I  =  62.1°,  €*95  =  5.0°  for  N  =  16, 
having  been  corrected  from  a  bedding  attitude  of  strike  =  050°  dip  =20°  to  original  dip 
of  strike  =  050°  dip  =10°  (strike  given  as  azimuth  clockwise  from  north  such  that  the  dip 
direction  is  to  the  right;  right-hand  rule).  This  indicates  21.6°±8.6°  clockwise  vertical- 
axis  rotation  (method  of  Demarest,  1983)  of  the  fault  block  on  which  these  sediments  lie. 
We  expect  these  sediments  to  average  secular  variation  because  the  thickness  spanned 
(12m)  contains  evidence  of  many  separate  depositional  events  each  being  a  few  tens  of 
centimeters  thick.  As  previously  stated,  the  fanglomerate  is  separated  from  the 
overlying  pyroclastics  and  dadte  by  an  angular  imconformity.  Given  that  both  parts  of 
the  section  are  of  similar  age,  we  lump  them  together  in  determination  of  locality  mean 
direction.  Separation  of  the  two  parts  of  the  section  results  in  an  even  larger  vertical 
axis  rotation  for  the  fanglomerate. 

3.4.2  Lava  Localities 

We  sampled  22  localities  of  basalts,  andesites,  and  dacites  of  the  Wild  Horse  Mesa 
spanning  ages  of  1.4  -  3.6  Ma  (see  Geochronology  section).  Of  this  group  several  were 
composed  of  many  separate  lava  flows  with  up  to  21  flows  represented.  The  hallmark 
of  these  thick  sections  however,  is  the  repetition  of  remanence  directions  and 
geochemistry  for  many  flows  within  stratigraphically  bounded  zones.  These  have  been 
interpreted  to  be  single  eruptive  units,  emplaced  in  a  short  span  of  decades  or  less  time. 
As  many  as  19  flows  are  observed  to  compose  a  single  eruptive  unit  (HS  locality.  Figure 
3.24),  usually  where  the  flows  pile  up  at  a  break  in  paleo-slope,  witiun  paleovalleys,  or 
near  the  eruptive  source  for  vents  located  on  low  paleoslopes.  This  results  in  only  one 
or  a  few  remanence  directions  being  represented  at  any  given  locality.  Therefore, 
averaging  secular  variation  is  very  difficult  for  the  lavas  of  this  region.  To  suffidently 
average  secular  variation,  a  combination  of  data  from  many  localities  would  be 
necessary,  and  the  resulting  loss  of  spatial  resolution  might  mean  liunping  rotated  with 
unrotated  blocks.  For  these  reasons  comparison  of  remanence  directions  of  individual 
eruptive  units  between  localities  yields  more  information  about  the  tectonics  of  the 
region. 

Data  in  Figure  3.24  and  3.25  illustrate  the  previous  points.  Depicted  are  flow  mean 
ChRMs  on  one  stereonet  for  each  sampling  locality  or  eruptive  umt  For  example,  the 
HS  locality  spans  nineteen  lava  flows  but  ttie  mean  flow  directions  of  these  nineteen  are 
statistically  indistinguishable,  indicating  eruption  in  a  period  shorter  than  that 
necessary  for  significant  secular  variation,  about  two  centuries.  Even  at  PC,  the  locality 
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representing  the  most  eruptive  units  (four),  this  is  insufficient  to  average  secular 
variation.  Only  the  AL  locality,  exhibiting  seven  distinct  remanence  directions  within 
the  Coso  Basin  andesite  unit,  may  sufficiently  average  secular  variation. 

Without  the  possibility  of  averaging  secular  variation,  we  use  the  spatial  extent  of 
eruptive  omits  and  the  high  precision  with  which  their  remanence  direction  can  be 
measured  to  detect  differences  in  remanence.  Figoire  3.26  shows  comparisons  between 
localities  of  unit  mean  ChRMs  of  eruptive  omits  shared  in  common.  Figure  3.26a  depicts 
mean  ChRM  by  locality  of  Petroglyph  Canyon  basalt  and  Black  Canyon  flows  for 
localities  BC,  PC,  LDC,  UDC,  FS.  In  the  geographic  reference  frame  (without  tilt 
correction)  locality  ChRMs  for  new  Basalt  of  Petroglyph  Canyon  are  mostly  statistically 
indistingoiishable  at  95%  confidence  while  remanence  directions  from  Black  Canyon 
flows  at  two  of  these  localities  are  indistinguishable  when  confidence  levels  are  raised 
to  98%.  On  the  other  hand  with  tilt  corrections  applied  (the  largest  being  only  10°), 
within-unit  directions  diverge,  indicating  that  at  least  part  of  the  tilt  of  lava  flow  tops 
results  from  original  dips  at  the  time  of  emplacement.  Therefore,  if  we  retain  the 
directions  with  no  tilt  correction  and  assume  that  some  of  these  flow-top  attitudes  stem 
from  original  emplacement  on  slopes,  there  is  no  relative  vertical  axis  rotation  of  these  5 
localities  relative  to  one  another.  Thus  if  one  locality  has  experienced  rotation,  all  must 
have. 

Figure  3,26b  shows  data  from  the  Cole  Flat  basalt  found  at  the  CF  and  FS  localities. 
In  geographic  and  tilt-corrected  coordinates  there  is  a  declination  difference  between 
locality  mean  directions.  In  geographic  coordinates  FS  is  rotated  11.6°  ±  6.9°  (method  of 
Demarest,  1983)  clockwise  about  a  vertical  axis  relative  to  CF,  while  tilt-corrected  this 
rotation  is  12,0°  ±  7.0°  (with  a  misfit  of  inclination  as  well).  As  will  be  shown  later,  this 
suggestion  that  FS  has  rotated  with  respect  to  a  "fixed"  CF  makes  most  sense. 
Remanence  data  from  Silver  Peak  basalt,  which  crops  out  at  SP  and  PC  localities,  are 
depicted  in  Figure  3.26c.  PC  is  clockwise  rotated  about  a  vertical  axis  relative  to  SP 
12,3°  ±  6.2°  (geographic)  or  13.0°  ±  6.2°  (tilt-corrected). 

To  summarize,  using  Cole  Flat  basalt  we  measmre  11.6°  clockwise  vertical-axis 
rotation  of  FS  relative  to  CF  locality.  Using  Petroglyph  Canyon  basalt,  FS  is  unrotated 
relative  to  PC  (as  well  as  relative  to  BC,  UDC,  and  LDC).  Using  Silver  Peak  basalt,  PC  is 
12.3°  clockwise  rotated,  about  a  vertical  axis,  relative  to  SP  locality.  Thus,  we  have  two 
corroborative  measures,  in  sense  and  magnitude,  of  clockwise  relative  rotation  of  five 
Wild  Horse  Mesa  localities  relative  to  two  at  its  margins. 

The  CHS  locality,  on  a  moxmtaintop  just  southwest  of  Coso  Hot  Springs,  is 
composed  of  poorly-preserved  outcrop  of  Black  Canyon  flows.  This  locality  is  on  the 
west  side  of  the  Coso  Wash  fault  zone  and  thus  perhaps  in  a  separate  tectonic  domain 
from  the  Wild  Horse  Mesa  and  Coso  Wash.  Further  west  there  is  abundant  evidence  of 
down  to  the  west  tilting  and  Whitmarsh  et  al.,  (1996)  found  no  net  vertical-axis  rotation 
since  the  Cretaceous.  Figure  3.26d  depicts  mean  unit  ChRMs  from  all  localities 
containing  Black  Canyon  flows.  There  is  no  available  tilt  correction  for  CHS  because  of 
poor  preservation.  However,  it  is  clear  from  the  difference  in  mean  locality  ChRM  for 
CHS  that  some  tilt  correction  is  required.  Therefore,  we  use  the  directional  difference 
between  CHS  and  Black  Canyon  flows  at  BC  to  determine  the  tOt  of  the  fault  block  on 
which  the  former  locality  sits.  Given  the  ChRM  for  CHS  presented  here  it  makes  little 
difference  whether  CHS  and  BC  have  both  rotated  about  a  vertical  axis,  or  only  BC  has. 
In  eitiier  case  the  tilt  required  to  restore  directional  correspondence  is  about  3()°±5° 
about  an  axis  with  azimuth  of  132°  and  a  plunge  of  between  15°  and  33°.  The  presence 
of  tilt  at  this  locality  is  the  significant  thing,  whereas  the  amount  of  tilt,  tilt  axis,  and 
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discerning  whether  BC  rotated  about  a  vertical  axis  independently  of  CHS  is  uncertain 
due  to  the  poor  preservation  of  outcrop  at  CHS.  Unlike  all  other  lava  localities 
described  in  this  report,  CHS  specimen  ChRMs  exhibited  a  large  dispersion.  The  small 
(X95  is  possible  for  the  site  because  of  a  large  number  of  samples  collected  (N=35),  with 
rejection  of  a  few  samples  exhibiting  ChRMs  suggesting  out-of-place  blocks  on  the 
outcrop  scale.  Therefore,  accuracy  of  the  precise  locality  direction  is  compromised.  In 
short  CHS  is  clearly  tilted  down  to  the  west,  by  anywhere  from  30°  to  50°  (the  latter 
magnitude  if  a  few  samples  are  rejected) 

In  the  area  around  Haiwee  Spring,  we  have  sampled  three  localities  of  Haiwee 
Spring  basalt  (localities  HS,  HSB,  and  HSD).  Figure  3.26e  presents  the  data,  but  neither 
the  geographic  nor  simplistic  tilt-corrected  coordinates  result  in  correspondence  in 
directions  for  tiiese  localities.  For  this  to  be  evidence  of  relative  vertical-axis  rotation 
they  must  at  least  exhibit  similar  inclinations.  This  is  only  the  case  for  HS  and  HSD  in 
the  imcorrected  case. 

There  are  several  possibile  explanations  for  this.  It  is  possible  that  although 
geochemically  the  same,  these  flows  may  record  different  times  in  Earth's  secular 
variation  in  magnetic  field.  This  is  improbable  because  of  the  thickness  of  the  unit 
(nineteen  flows  at  HS),  close  proximity  of  the  localities  and  the  absence  of  remanence 
variation  within  a  locality.  These  all  suggest  rapid  voluminous  eruption  that 
simultaneously  covered  the  small  vicinity  of  these  localities.  It  is  possible  that 
unaccoimted-for  original  dips  cause  part  of  the  mismatch  with  the  additional  possibility 
of  relative  rotation  of  localities.  For  example,  simplistic  tilt-correction  of  HS  and 
correction  of  HSB  to  an  original  dip  of  about  10°  southwards  results  in  equivalent 
inclinations  and  a  clockwise  rotation  of  HSB  relative  to  HS.  The  number  of  degrees  of 
freedom  on  such  solutiojis  is  infinite  -  one  must  solve  for  original  dip  of  each  pair  of 
localities  plus  relative  vertical-axis  rotation  between  them  given  the  hard  constraint  of 
the  remanence  directions  and  the  guideline  of  flow  "attitudes".  More  information 
about  the  flow  direction  from  knowledge  of  the  eruptive  source  would  help  in 
constraining  the  solution.  In  short  Haiwee  Spring  localities  suggest  original  dips  and 
perhaps  relative  vertical-axis  rotation,  but  a  quantitative  measure  is  premature  pending 
additional  data. 

Coso  Basin  andesites  first  recognized  by  Whitmarsh  (1996)  to  be  Pliocene  (map  unit 
Tbb  "basalt"  of  Coso  Basin,  3.58±0.86  Ma)  in  age  are  the  only  lavas  in  the  field  area 
having  any  possibility  of  yielding  an  absolute  rotation  through  averaging  of  secular 
variation.  Along  the  Coso  Wash  fault  zone,  the  AL  locality,  1.5  km  north  of  Airport 
Lake's  western  end,  contains  5  flows,  a  dike,  and  a  sill  ( Figiure  3.26f)  each  recording  a 
different  remanence  direction.  The  flows  have  been  tilted  to  the  west  back  toward  the 
presumed  eruptive  source,  a  near-surface  intrusive  body  of  andesite.  When  tilt- 
corrected  to  )tield  an  assiuned  original  dip  of  10°  to  the  southeast  the  locality  5tields  a 
mean  remanence  direction  indistinguishable  from  geographic  south  though  somewhat 
shallow.  The  shallow  inclination  indicates  that  either  tiie  tilt  correction  is  inappropriate 
or  that  secular  variation  was  not  fully  averaged  out. 

Several  Plio-pleistocene  lavas  indicate  no  relative  rotation  between  the  localities 
where  they  occur.  These  negative  results  constrain  the  spatial  pattern  or  timing  of 
vertical-axis  rotation  and  tilting  in  the  study  area.  Duffield  and  Bacon  (1981)  mapped 
Basalt  of  Renegade  Canyon  and  subdivided  it  into  multiple  subunits,  QTbrj-QTbrg.  We 
have  sampled  QTbrj  and  QTbr4  at  two  localities  each  and  found  no  statistically 
distinguishable  directional  difference  within  each  pair  of  directions  (Figures  3.26g,  h). 
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Either  differential  rotation  between  these  localities  has  been  insufficient  during  the 
period  measured  for  each  subunit  or  no  differential  rotation  has  occurred  in  this  region. 
Note  that  QTbrj  at  the  Renegade  Canyon  locality  exhibits  two  distinct  directions  for 
different  distal  lava  flows.  This  suggests  either  a  sampling  error,  that  secular  variation 
was  very  fast  during  their  emplacement  or  that  long  hiatuses  occurred  between 
emplacement  of  the  two  flows. 

We  sampled  Basalt  of  Coso  Wash  (localities  CMA  and  CMC  from  Qbc,  Duffield  and 
Bacon,  1981),  in  order  to  place  constraints  on  the  differential  tilting  rate  across  the  Coso 
Wash  fault  zone.  Four  kilometers  north  of  this  unit  is  the  Coso  Hot  Spring/Black 
Canyon  locality  pair  with  its  30°-50°  differential  tilt  since  3.4  Ma.  Two  kilometers  south 
is  the  Airport  Lake  andesite  which  has  tilted  10°-20°  since  3.6  Ma.  We  sampled  Qbc 
(dated  at  1.1  Ma  by  Duffield  et  al,  1980,  ca.  0.97  Ma  according  to  our  geochronology) 
with  the  intent  of  measuring  some  intermediate  amoimt  of  relative  tilt  across  the  Coso 
Wash  fault  zone.  We  found  no  tilt  (or  vertical-axis  rotation)  within  an  imcertainty  of 
only  a  few  degrees.  This  constrains  initiation  of  tilting  on  the  Coso  Wash  to  have 
predated  1  Ma.  Otherwise  a  relative  tilt  would  probably  be  observed  between  localities 
CMA  and  CMC 

3.5  Geochronology 

A  summary  of  oxu:  geochronologic  data  is  presented  in  Table  3.1.  Two  of  the 
analyses  provided  poor  results  because  of  the  character  of  the  lavas.  The  remainder  of 
lavas  sampled  pelded  satisfactory  and  sometimes  imexpected  results.  For  example,  a 
sample  from  RC  locality  which  correlates  paleomagentically  and  geochemically  with 
QTbrj  was  determined  to  be  1.37±0.04  Ma,  while  Duffield  et  al  (1980)  assigned  an  age  of 
1.75±0.1  Ma  based  on  K/Ar  dating  on  samples  from  locality  RCE.  In  addition  *Ar/^’Ar 
geochronology  from  the  base  of  RC  locality  yielded  an  age  of  about  2.6  Ma, 
substantially  older  than  the  other  radiometric  constraints  for  rocks  around  Renegade 
Canyon.  Remaining  geochronologic  determinations  are  of  good  quality  and  yielded 
ages  similar  to  those  expected. 

The  most  significant  result  of  our  new  geochronologic  data  is  that  they  suggest  that 
eruption  of  Wild  Horse  Mesa  lavas  spanned  a  shorter  period  than  previously  thought. 
Duffield  et  al.  (1980)  considered  these  eruptions  to  have  spanned  3-4  Ma,  though  they 
acknowledge  most  eruptions  postdated  3.7  Ma.  Some  of  their  K/  Ar  dating  included 
ages  older  tt\an  4  Ma  (e.g.  their  localities  numbered  38  and  47).  Our  combined 
geochronology  and  stratigraphic  constraints  show  that  eruptions  in  the  Wild  Horse 
Mesa  began  arotmd  3.6  Ma  and  were  largely  completed  by  3.0  Ma  (Figure  3.27),  when 
volcanic  activity  shifted  to  the  Haiwee  Reservoir  area.  Whereas  Wild  Horse  Mesa 
volcanics  were  voluminous  and  of  a  "short"  duration,  volcanic  activity  in  the  Renegade 
Canyon  area  was  of  smaller  volume  and  much  more  distributed  in  time,  spanning  2.6  to 
1.4  Ma.  The  short-dxiraition,  voluminous  Wild  Horse  Mesa  volcanics  suggest 
relationship  with  regional  extensive  volcanism  from  aroxmd  the  same  period,  now 
being  attributed  to  foimdering  of  an  ecologite  root  from  beneath  the  Sierra  Nevada 
(Manley  et  al.,  2000). 

4.0  Discussion 

The  Pliocene  lavas  and  sediments  sampled  in  Wild  Horse  Mesa  and  Coso  Wash  all 
indicate  clockwise  vertical-axis  rotation  of  the  fault-bounded  blocks  sampled.  The  Coso 
Wash  fanglomerate  locality  reveals  absolute  rotation  relative  to  paleo-geographic  north. 
The  five  Wild  Horse  Mesa  lava  localities  (BC,  PC,  LDC,  UDC,  and  FS)  are  little  changed 
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relative  to  one  another  but  clockwise  rotated  relative  to  CF  and  SP.  The  uncertainties  on 
the  measurements  presented  here  are  not  small  enough  to  discern  whether  there  is 
spatial  variability  in  clockwise  rotation  within  the  clockwise-rotated  region. 

Figure  4.1  shows  a  simplified  map  of  the  geometry  of  rotated  blocks  relative  to 
boundaries  of  the  rotating  region.  The  geometry  of  the  rotating  region  is  more  complex 
than  models  for  block  rotation  usually  presented  in  the  scientific  literature  (e.g., 
Luyendyk,  1991;  Nur  et  al.,  1993,  Unruh  et  al.  1996),  because  the  boundaries  of  the 
rotating  zone  are  not  parallel.  These  boimdaries  are  the  Coso  Wash  (normal)  fault  zone 
and  the  herein  defined  Wild  Horse  Mesa  fault  zone.  The  geometry  of  this  latter 
structure  is  estimated  using  (1)  the  geometry  of  mapped  fault  traces,  (2)  the  edge  of 
mapped  Mesozoic  basement  of  the  Cole  Flat  area,  and  (3)  the  fact  that  tius  feature 
separates  a  domain  of  vertical-axis  rotation  from  a  larger  domain  that  has  not 
experienced  similar  rotation  to  the  east  (Schweig,  1989;  Ron  and  Nur,  1996).  As  defined 
above  it  appears  as  a  right-stepping,  northwest-striking  string  of  faults  nearly  optimally 
oriented  to  accomodate  regional  dextral  shear.  In  fact,  the  M  5.2,  March  1998 
earthquake  occurred  sub-parallel  to  and  nearby  this  structiure. 

We  hypothesize  that  block  rotation  at  the  near-surface  at  Wild  Horse  Mesa  and  Coso 
Wash  accommodates  regional  dextral  shearing.  Block  long-axes  are  oriented  between 
due  north  and  N20°E  while  geodeticaUy-determined  dextral  motion  across  the  Eastern 
California  shear  zone  at  this  latitude  (Airport  Lake  fault  zone)  is  oriented  at  about 
N40°W  (McClusky  et  al.,  2001).  Similarly,  the  geodeticaUy-determined  azimuth  of 
motion  of  the  Sierra  Nevada  microplate  relative  to  stable  North  America  is  N43°W  at 
this  latitude  (Dixon  et  al.  2000).  Consequently,  fault-botmded  blocks  of  the  Wild  Horse 
Mesa  are  oriented  more  favorably  for  accomodating  the  regional  shear  by  rotation  than 
either  pure  dextral  slip  or  pure  dip-sUp  on  the  faults  between  them  (e.g.,  Dewey  2002). 

The  orientation  of  die  WUd  Horse  Mesa  and  Coso  Wash  fault  zones  and  the  fact  that 
these  boimdaries  are  not  parallel  requires  net  thinning  in  between  as  part  of  the 
regional  dextral  shearing.  Early  pinned  block  models  of  rotation  required  opening 
perpendicular  to  the  parallel  basin  boimding  faults  during  transtension  (e.g.  Luyendyk, 
1991).  At  WUd  Horse  Mesa  opening  of  the  basin  is  accomplished  passively  by  the  basin 
geometry.  Furthermore,  in  pinned  block  models  finite  rotation  requires  obUque  slip 
with  either  a  component  of  normal  or  reverse  faulting  between  blocks,  depending  upon 
the  geometry.  In  the  WUd  Horse  Mesa  only  dip-slip  offsets  are  apparent  because  of  the 
lack  of  horizontal  offset  markers.  This  crustal  thini^g,  as  demonstrated  by  inverse 
modeling  of  seismogenic  strain  at  shaUow  crustal  levels  (Figure  3.1),  is  therefore 
compatiWe  with  vertical  axis  rotation.  Moreover,  the  outcrop-scale  brittle  faults 
exposed  on  WUd  Horse  Mesa  record  components  of  crustal  thinning  and  horizontal 
shearing,  both  of  which  are  characterized  by  subhorizontal,  west  southwest-east 
northeast  maximum  stretching  (dj).  These  strain  geometries  are  consistent  with  the 
seismogenic  strain  geometries  as  weU  as  expectations  of  transtension  in  the  regional- 
scale  right  releasing  step  in  the  dextral  faults  that  mark  the  eastern  California  shear 
zone  in  this  area  (e.g..  Figure  2.1,  and  Unruh  et  al.  2002).  Therefore,  these  data  present  a 
unified  picture  of  khiematic  behavior  of  the  WUd  Horse  Mesa. 

Block  models  are  useful  simpUfications  because  tiiey  permit  calculation  of  the 
minimum  total  dextral  sUp  necessary  to  produce  the  observed  rotations.  Across  such  a 
region,  more  direct  measurements  of  displacement  are  often  impossible.  Figure  4.2 
depicts  such  a  calculation  of  minimum  dextral  sUp  on  the  basin-boimding  faults  for  the 
observed  rotations,  asstuning  blocks  pinned  to  the  boundary  faults.  Typically,  such 
calculations  require  block  length,  width  and  orientation  relative  to  strike-slip  boundary 
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fatilts.  Assuming  tiiat  the  blocks  are  much  longer  than  they  are  wide,  which  is 
consistent  with  ihe  mapped  fault  traces,  the  exact  width  of  the  blocks  becomes  less 
important  to  this  calculation.  This  is  especially  so  in  the  Wild  Horse  Mesa  case  because 
normal  faulting  results  in  apparent  widening  of  blocks  in  map  view  as  block  footwalls 
are  exhumed  across  inclined  faults.  In  addition,  we  allow  the  blocks  to  deform 
internally  by  slip  on  faults  of  the  type  examined  in  this  study  (i.e.,  we  treat  ihe  shallow 
crust  as  a  micropolar  continuum).  This  view  is  supported  by  the  fact  that  closely- 
spaced  outcrop-scale  faults  have  accommodated  crustal  thiiming,  horizontal  shearing 
and  crustal  thickening.  Using  an  estimated  block  length  of  10  km  obtained  by  using 
crust  bounded  by  only  the  largest  normal  faults  in  the  study  area,  we  calculate  a 
nunimum  dextral  offset  of  1.5  km  along  the  Wild  Horse  Mesa  basin-boimding  fault 
zone.  The  calculated  minimum  offset  is  easily  achievable  given  the  current 
geodetically-determined  slip  rates  across  the  basin  of  5.3  mm/yr  (the  rate  of  Airport 
Lake  fault  zone,  McClusky  et  al.,  2001).  The  total  minimiun  offset  is  accotmted  for  in 
only  a  few  htmdred  thousand  years,  given  this  as  a  steady  rate,  despite  current 
imderstanding  that  dextral  shear  initiated  in  the  region  as  recently  as  ca.  2  Ma 
(Monastero  et  al.  2002).  This  implies  additional  slip  not  accoimted  for  in  our  calciilation 
arising  either  from  impinning  the  rotated  blocks  or  from  other  structures  outside  the 
Wild  Horse  Mesa. 

It  is  important  to  understand  ttiat  rotating  blocks  of  Wild  Horse  Mesa  are  an  upper- 
crustal  feature  probably  riding  above  one  or  a  few  dextral  master  faults  at  depth,  as 
indicated  by  the  differing  character  of  seismogenic  strain  with  depth.  In  addition, 
significant  complexity  is  not  accounted  for  in  the  model,  such  as  tire  possibility  of 
curved  faults,  varying  tilt  magnitude  on  different  blocks  and  block  rotation  in  the  area 
west  of  Coso  Wash  fault  zone.  Nonetheless,  our  block  model  is  instructive  and  retains 
the  most  important  features  of  the  study  area. 

Within  the  Wild  Horse  Mesa,  field  observations  indicate  that  some  fault-bounded 
blocks  have  tilted.  In  particular  some  lava  flow  tops,  which  parallel  the  ground  surface 
across  most  of  the  area,  display  substantial  topographic  slopes.  Adjacent  blocks  tilt  in 
different  directions  and  back  toward  eruptive  sources  from  which  lavas  emanated  in 
some  cases.  Our  sampling  localities  have  inadvertantly  avoided  such  blocks  for  the 
most  part.  The  little  tilt  present  at  PC,  BC,  UDC,  and  FS  can  be  attributed  to  original  dip 
of  the  lavas  as  they  were  emplaced  on  a  paleoslope,  as  supported  by  the 
correspondence  of  remanence  directions  for  these  localities  despite  the  apparent  tilt. 
This  is  a  common  phenomenon  in  other  studies.  In  the  case  of  the  three  Haiwee  Springs 
localities  (HS,  HSB,  HSD),  correction  for  large  tilts  results  in  great  discrepancies  in 
remanence  direction  within  a  single  eruptive  unit.  This  is  probably  best  attributed  to 
steep  original  dips  coupled  with  differential  rotation  between  localities.  We  suggest 
tiiat  ihe  tilting  of  fault  blocks  in  the  Wild  Horse  Mesa  is  not  a  significant  feature.  For 
example  ihe  tilts  in  opposite  directions  (some  east,  some  west)  in  southern  Wild  Horse 
Mesa  suggest  a  net  canceling  out  of  at  least  some  this  horizontal  axis  rotation  in  that 
region.  Our  paleomagnetic  sampling  of  large  stable  blocks  within  the  mesa  indicates 
little  net  tilt.  Sandbox  modeling  (Dooley  and  McClay,  1997)  of  step  overs  in  strike  slip 
fault  systems  supports  this  and  indicates  that  net  vertical-axis  rotation  is  the  more 
significant  and  regioneJ  feature. 

However  tilt  on  the  west  side  of  the  Coso  Wash  fault  zone  appears  to  be  generally 
down  to  the  northwest.  Coso  formation  between  Haiwee  Reservoir  and  the  geothermal 
field  is  mostly  tilted  in  this  sense.  Paleomagnetic  data  from  CHS  and  lava  flow 
attitudes  from  AL  and  HSD  indicate  large  tilts  of  20°  to  35°  or  more  along  the  west  side 
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of  the  Coso  Wash  fault  zone.  Whitmarsh  et  al.  (1996)  suggested  similar  tilting  of  30°  or 
more  near  Rose  Valley  since  the  Cretaceous.  Altogether  Ihis  indicates  widespread 
tilting  down  to  the  northwest  of  Mesozoic  basement  blocks  between  the  Coso  Wash 
fault  zone  and  Rose  Valley /Haiwee  Reservoir.  Whether  vertical  axis  rotation  has 
occurred  in  this  region  is  still  debatable  though  Whitmarsh's  impublished  data  are  to 
the  contrary. 

5.0  Conclusion 

The  Wild  Horse  Mesa  exhibits  12.0°±4.6°  to  21.6°±10.7°  of  clockwise  vertical-axis 
rotation  since  about  3  Ma  of  fault  boimded  blocks  as  demonstrated  paleomagnetically. 
Abundant,  probably  oblique  sinistral-normal,  north-northeast  striking  faults 
accommodate  block  rotation  near  the  surface  between  the  left-stepping  Coso  Wash 
normal  fault  zone  and  the  newly-defined,  right-stepping,  dextral  Wild  Horse  Mesa  fault 
zone.  Clockwise  block  rotation  is  consistent  with  our  results  of  inverse  modeling  of 
fault  slip  data.  Moreover,  statistically-significant  relative  vorticity  is  suggested  by  our 
inverse  modeling  of  eartiiquakes  from  5-8  km  at  Wild  Horse  Mesa  (Figure  3.2).  Given 
the  predominance  of  regional  right  slip,  and  the  fact  that  the  earthquake  hypocenters 
define  a  northwest  striking  seismic  source  zone,  the  sense  of  block  rotation  for  these 
seismic  events  is  interpreted  to  be  clockwise.  Absolute  vertical  axis  rotations  of  12° 
yield  a  minimum  dextral  offset  across  the  Wild  Horse  Mesa  of  1.5  km  along  a  strike  of 
N30°W,  subparallel  to  the  Wild  Horse  Mesa  fault  zone. 

The  dextral  fault(s)  at  depth  merge  with  or  step  over  to  the  Airport  Lake  fault  zone 
to  the  south.  We  suggest  that  to  the  north  this  strike-slip  fault  steps  across  Upper 
Centennial  Flat  to  Lower  Centennial  Flat  and  the  east  side  of  Owens  Lake  where 
substuface  faults  have  been  mapped  by  seismic  reflection  (Neponset  Geophysical  and 
Aquila  Geosciences,  1997).  This  semi-continuous  structure  is  one  strand  of  the  Eastern 
California  shear  zone  lining  the  Owens  Valley  fault  zone  with  the  Airport  Lake  fault 
zone,  and  ultimately  with  structures  south  of  the  Garlock  fault.  The  Mesozoic  basement 
high  to  the  west  of  this  localized  strand  of  the  Eastern  California  shear  zone  is  tilted  as 
many  separate  blocks  to  the  northwest.  This  tilting  could  be  the  result  of  upper  plate 
disruption  above  a  master  Sierra  Nevada  frontal  or  other  fault  zone  (e.g.,  Monastero 
1997).  Alternatively,  a  through-going  structure  may  not  be  necessary  or  even  possible 
beneath  the  geothermal  field  since  the  magma  chamber  there  probably  accommodates 
extension. 

In  summary  we  have  shown  that  careful  mapping  and  geochemical  fingerprinting  of 
lava  flows,  in  conjunction  with  paleomagnetic  analyses  and  radiometric  dating,  provide 
important  constraints  on  volcano-stratigraphy.  These  constraints  have  been  integrated 
with  analyses  of  contemporary  non-recoverable  strain  (i.e.,  earthquakes)  and  time- 
integrated  non-recoverable  strain  (i.e.,  exposed  faults),  to  derive  a  better  imderstanding 
of  the  kinematic  evolution  of  the  Wild  Horse  Mesa  and  adjoining  Coso  geothermal  field. 
The  result  is  (1)  a  conceptual  model  for  the  geometry  of  some  of  the  major  structures 
that  accommodate  transtension  in  the  region,  and  (2)  an  estimate  of  the  magnitude  of 
slip  required  to  accotmt  for  the  observed  finite  rotations  of  crustal  blocks,  the  geometry 
of  contemporary  strain  and  the  geometry  of  long-term  fault-related  strain. 
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Table  1.1  Activity  Summary  By  Task 


Task 

Description 

Activity  dxiring  quarter  (yes/no) 

1 

Inversion/evaluation  of  earthquake  data 

No 

2 

Micropolar  inversion  of  field-measured  fault  data 

Yes 

3 

Radiometric/Paleomagnetic  analyses 

Yes 

4 

Integration  and  interpretation 

Yes 

Purchase  Request  Number:  N68936-01-C-0094 
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Table  2.1  Field-measured  shear  plane/slip  data  from  Wild  Horse  Mesa. 


Purchase  Request  Number:  N68936-01-C-0094 


25 


Purchase  Pjequest  Number:  N68936-01-C-0094 


26 


u 

232 

’  282.6 

70.9 

2 

213 

u 

35 

■bbee 

138.5 

72.5 

8 

213 

u 

19 

30 

118.2 

29.7 

8 

213 

u 

6 

88 

156 

86 

8 

213 

u 

32 

67 

o 

CO 

8 

213 

u 

50 

66 

178.6 

60.3 

8 

213 

u 

255 

69 

336.7 

68.8 

8 

213 

u 

45 

77 

49 

16 

4 

u 

300 

74 

12.39 

73.26 

2 

1 

u 

2 

80 

34 

72 

2 

u 

149 

43 

326 

4 

2 

u 

2 

82 

48 

79 

2 
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^  "a"  and  "b"  quality  not  differentiated,  see  text  for  discussion  of  use  of  "c"  quality  data 
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Table  3.1  "Ar/^’Ar  Geochrotiology. 


Erupfave 

Unit 


Haiwee 

Springs 

basalt 


Location 


QTbri  (?) 


QTbri  (?) 


QTbij 


Black 

Canyon 

flows 


Basalt  of 
Coso 
Wash 


Petroglyph 

Canyon 

basalt 


Petroglyph 

Canyon 

basalt 


Petrogl3^h 

Canyon 

basalt 


Petroglyph 

Canyon 

basalt 


Petroglyph 

Canyon 

basalt 


Petroglyph 

Canyon 

basalt 


Cole  Flat 
basalt 


CWA 


CMA 


Sample 

Name, 

Flow# 


HS19, 


Method 

and 

Material 


Isochron, 


flow  19*  Groudmass 


Tdo, 
Dnffield 
&  Bacon 
unit  Tdo 


RC14, 
flow  14t 


RC14, 
flow  14t 


RCOl, 
flow  It 


BC18, 
flow  18 


Isochron, 

Groudmass 


Plateau, 

Groudmass 


Comments  Age  (Ma) 
and  2  a 
error 


Argon  recoil,  >2.96±0.03 
age  is  a 
minumum 


Excess  <2.83±0.02 
Argon,  age  is 
a  maximum 


2.62±0.03 


Plateau, 

Groudmass 


Plateau, 

Groudmass 


2.63±0.04  I  2.62±0.02 


1.37±0.04 


BD0201A,  Isochron, 
flow  2  J  Whole 
Rock 


BD0301A,  Isochron, 
flow  3 1  Whole 
Rock 


BD0303A, 
flow  3 1 


3.41±0.04 


Previously  c.a.  0.97 
mis-dated  or 
mis-identified 


3.19±0.04 


3.20±0.04 


BD1403A, 
flow  14  ^ 


BD0802A, 
flow  8^ 


3.20±0.03  I  3.22±0.02 


3.25±0.03 


3.25±0.05 


3.26±0.04 


3.51±0.04 


3.50±0.03 
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Cole  Flat  I  FS  BD1204A,  Isochron,  3.49±0.05 

basalt  flow  12  $  Whole 

_ _ Rock _ 

*  numbered  from  bottom  of  stratigraphic  section,  t  numbered  from  top  of  stratigraphic 
section,  f  numbered  irregularly 
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Figure  2.1.  Seismically-defined  right-releasing  step  in  the  boundary  between  the  Sierran  microplate 
and  the  southern  Basin  and  Range  Province  in  the  area  of  the  Wild  Horse  Mesa  and  Coso  Geothermal 
Field  (Indicated  by  the  six-sided  figure).  The  rectangle  approximately  outlines  the  seismic  events 
described  in  Figures  3.1  and  3.2.  The  trend  of  subhorizontal  di  and  ds  for  best-fitting  models  for 
events  shown  in  light  grey  are  indicated  by  heavy  and  light  lines,  respectively.  The  best-fitting  great 
circles  for  these  event  epicenters  are  subparallel  to  the  trajectory  of  the  Sierran  microplate  with  respect 
to  North  America  (Dixon  et  al.,  2000)  shown  with  the  large  arrow. 
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\  \  Owens 

I  \  \  Lake 
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Figure  2.2  >  Geologic  and  Sampling 
Locality  Map.  Tertiary  and  Mesozoic 
basement  rocks  of  greater  Coso 
Range  are  shown.  Geology  is 

\  drawn  from  Duffield  and 

Bacon  (1981)  and  Stinson 
(1977a&b). 
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Figure  2.3  -  Geologic  map  of  the  Coso  Volcanic  field  showing  approximate  locations  where 
brittle  faults  were  examined  as  white  bands  within  the  Pliocene  volcanics  (solid  grey). 
Canyons  labeled  as  follows:  BCy  =  Black  Canyon,  DCy  =  Dead  End  Canyon,  Iwh  =  lower 
Wild  Horse  Mesa,  PCy  =  Petroglyph  Canyon  and  SCy  =  Sheep  Canyon,  uwh  =  upper  Wild 
Horse  Mesa. 
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Figure  2.4.  Example  of  shear-sense  indicators  preserved  on  outcrop-scale 
fault  surfaces.  Pencil  at  bottom  points  in  the  direction  of  the  missing  block. 
Lineated  subsidiary  shear  surfaces  (R  planes)  intersect  this  primary  slip 
surface  (Y  shear)  creating  two  lunate  features  above  and  to  the  right  of  the 
pencil  tip  (both  somewhat  shaded). 
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Figure  2.5.  Tangent  lineations  (left)  and  kinematic  P  and  T  axes  (right)  for  outcrop-scale  brittle  faults  at 
upper  and  lower  Wild  Horse  Mesa  (top  and  bottom,  respectively).  All  plots  are  lower-hemisphere,  equal- 
area  projections.  The  tangent  lineation  plots  show  poles  fault  planes  with  arrows  indicating  the  direction  of 
motion  of  the  footwall  blocks. 
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Wild  Horse  Mesa  Seismogenic  Strain 

Events  0-5  km  deep 


25  random  restarts 
2000  bootstrap  models 
D  =  0.596  (+0.114  -0.236) 
W  =  0.456  (+0.283  -0.232) 


best-fit  to  full  data  set: 
di  =  258.21  8.52 
d2  =  167.94  1.77 
d3  =  246.29  -81.30 
D  =  0.61 


Figure  3.1 .  (a)  Kinematic  P  and  T  axes  (with  separately  W  =  0.44 

determined  Kamb  contours)  for  earthquakes  from  0-5  km  in  V  =  -1 . 1 0 

depth  at  Wild  Horse  Mesa.  All  plots  are  lower-hemisphere,  mean  misfit  =  4.28 

equal-area  projections,  (b)  D  versus  W  plot,  (c)  Final  positions  D  not  statistically  different 
and  Kamb  contours  for  di  and  da  for  bootstrap  models  than  0.5 

(contoured  separately,  (d)  Strainball  plot  of  the  best-fitting 
strain  tensor  for  the  full  data  set. 
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Wild  Horse  Mesa  Seismicity  Strainball 

Events  5-8  km  deep 


0 

-1  0  1 

w 

Number  of  points  =  2000 
Bin  size  (D.  W]  =  [0.05, 0.05] 

Uniform  density  =  2.5  points,  or  0.125%,  per  bln 
Contours  at  zero,  uniform,  and  Intervals  of  39  times  uniform  density 
Peak  value  =  197  times  uniform  =  493  points  per  bln 


25  random  restarts 
2000  bootstrap  models 
D  =  0.569  (+0.204  -0.189) 
W  =  0.647  (+0.200  -0.426) 


best-fit  to  full  data  set: 
di  =  255.88  20.98 
d2  =  233.55  -67.48 
d3  =  342.87  -7.81 
D  =  0.57 


Figure  3.2.  (a)  Kinematic  P  and  T  axes  (with  separately  W  =  0.66 

determined  Kamb  contours)  for  earthquakes  from  5-8  km  in  V  =  0.19 

depth  at  Wild  Horse  Mesa.  All  plots  are  lower-hemisphere,  mean  misfit  =  4.94 

equal-area  projections,  (b)  D  versus  W  plot,  (c)  Final  positions  mean  misfit  =  4.28 

and  Kamb  contours  for  di  and  da  for  bootstrap  models  D  not  statistically  different 

(contoured  separately),  (d)  Strainball  plot  of  the  best-fitting 
strain  tensor  for  the  full  data  set. 
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25  random  restarts 
2000  bootstrap  models 
D  =  0.671  (+0.276  -0.340) 
W  = -0.011  (+0.407  -0.240) 


Figure  3.3.  (a)  Tangent  lineations  and  (b)  kinematic  P  and  T 
axes  for  outcrop-scale  brittle  faults  at  lower  Wild  Horse  Mesa 
that  display  a  shallow  west  northwest-plunging  P  axis 
maximum.  All  plots  are  lower-hemisphere,  equal-area 
projections,  (c)  Final  positions  and  Kamb  contours  for  di  and 
da  for  bootstrap  models  (contoured  separately),  (d)  Strainball 
plot  of  the  best-fitting  strain  tensor  for  the  full  data  set. 


best-fit  to  full  data  set: 
di=  171.71  37.16 
d2  =  234.97  -30.70 
da  =  297.62  37.74 
D  =  0.70 
W  =  -0.02 
V  =  -0.04 
mean  misfit  =  6.34 
D  and  W  not  statistically  different 
than  0.5  and  0.0,  respectively 
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25  random  restarts 
100  bootstrap  models 
D  =  0.935  (+0.065  -0.483) 
W  =  0.503  (+0.271  -1.091) 


Figure  3.4.  (a)  Tangent  lineations  and  (b)  kinematic  P  and  T 
axes  (with  separately  determined  Kamb  contours)  for  outcrop- 
scale  brittle  faults  at  upper  and  lower  Wild  Horse  Mesa  that 
accommodate  normal  faulting.  All  plots  are  lower-hemisphere, 
equal-area  projections,  (c)  Final  positions  and  Kamb  contours 
for  di  and  da  for  bootstrap  models  (contoured  separately),  (d) 
Strainball  plot  of  the  best-fitting  strain  tensor  for  the  full  data 
set. 


best-fit  to  full  data  set: 
di  =236.43  -8.16 
d2  =  148.23  12.38 
da  =  113.82  -75.10 
D  =  0.94 
W  =  0.50 
V  =  -1.66 
mean  misfit  =  5.27 
D  and  W  not  statistically  different 
than  0.5  and  0.0,  respectively 
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Crustal  Thickenina 


25  random  restarts  best-fit  to  full  data  set: 

100  bootstrap  models  - 119.94  -77.89 

D  =  0.073  (+0.695  -0.073)  d2  =  106.26  11.69 

W  =  0.191  (+0.640  -0.986)  d3  =  16.83  -2.77 

D  =  0.05 

_  W  =  0.25 

Figure  3.5.  (a)  Tangent  lineations  and  (b)  kinematic  P  and  T  V  =  0.94 

axes  (with  separately  determined  Kamb  contours)  for  outcrop-  mean  misfit  =  5.61 

scale  brittle  faults  at  upper  Wild  Horse  Mesa  that  accommodate  d  and  W  not  statistically  different 
thrusting.  All  plots  are  lower-hemisphere,  equal-area  than  0.5  and  0.0,  respectively 

projections,  (c)  Final  positions  and  Kamb  contours  for  di  and 
da  for  bootstrap  models  (contoured  separately),  (d)  Strainball 
plot  of  the  best-fitting  strain  tensor  for  the  full  data  set. 
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Upper  Wild  Horse  Mesa 
Horizontal  Shearing 


95''’o  Confidence  Limits  to?  D  and  W 


0 

w 


Number  of  points  =  1900 
Bin  size  [D,  W]  =  [0.1, 0.1] 

Uniform  density  =  9.5  points,  or  0.5%,  per  bin 
Contours  at  zero,  uniform,  and  intervals  of  10  times  uniform  density 
Peak  value  =  50  times  uniform  =  475  points  per  bln 


20  random  restarts 
2000  bootstrap  models 
0  =  0.152  (+0.615  -0.152) 
W  =  -0.550  (+0.895  -0.394) 


prolate  tensor  with 
non-zero  relative 
vorticity  suggested 
at  90%  confidence 


Figure  3.6.  Results  for  outcrop-scale  brittle  faults  at  upper  Wild  Horse  Mesa  that  accommodate 
horizontal  shearing.  Same  plots  as  Figures  3.4  and  3.5  except  that  plots  of  D  versus  W  are 
included  for  the  95%  and  90%  (in  box)  confidence  levels  (c).  See  text  for  discussion  of  statistics 
and  interpretation  of  statistics  at  90%  confidence  (plot  at  right). 
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Figure  3.7  -  Combined  Trace  Element  XRF  and  Paleomagnetic  Fingerprinting,  Pliocene  Wild 
Horse  Mesa  and  Drillhole  35-2.  Magnetic  polarity  is  indicated  by  solid  symbols  for  normal,  open 
symbols  for  reversed,  charcters  for  unknown  polarity.  The  use  of  paleomagenmtism  and  geochmis- 
try  allows  rigorous  separation  of  units,  as  indicated 
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Figure  3.8  -  Combined  Trace  Element  XRF  and  Paleomagnetic  Fingerprinting,  Plio-Pleistocene 
Lavas  of  Renegade  Canyon  Area.  Magnetic  polarity  is  indicated  by  solid  symbols  for  normal,  open 
symbols  for  reversed,  charcters  for  unknown  polarity.  The  use  of  paleomagnetism  and  geochemistry 
allows  rigorous  separation  of  units,  as  indicated.  RC  bottoms  flows  are  tentatively  correlated  with 
QTbrj  without  paleomagnetic  data. 
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Figure  3.9  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Silica  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.10  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Ti  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3,11  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  A1  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.12  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Fe  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.13  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Mn  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unlcnown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.14  •  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Ca  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.15  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Na  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.16  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  K  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  synrbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.17  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Phosphate  vs.  Mg.  Magnetic 
polarity  is  indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and 
characters  for  unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  sam¬ 
ples  deriving  from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.18  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Nb  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.19  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Zr  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35~2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.20  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Y  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.21  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Sr  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.22  -  Combined  XRF  and  Paleomagnetic  Fingerprinting,  Rb  vs.  Mg.  Magnetic  polarity  is 
indicated  by  solid  symbols  for  normal,  open  or  lightly  filled  symbols  for  reversed,  and  characters  for 
unknown  polarity.  These  data  corroborate  eruptive  unit  correlations  except  for  samples  deriving 
from  drillhole  35-2  which  probably  have  been  hydrothermally  altered. 
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Figure  3.23  -  Mean  Directions  of 
Sedimentary  Localities  shown 
with  blue  stars.  A.  WHB  locality 
is  not  rotated  within  the  rather 
large  0t95.  Thus  additional  work 
may  yet  reveal  significant 
rotation.  B.  Two  possible  tilt 
corrections  necessary  to  restore 
CWA  fanglomerate  to  its  origi¬ 
nal  dip  are  depicted.  In  either 
case  a  clockwise  vertical-axis 
rotation  is  revealed. 
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tiii-correctedl  Figure  3.24  -  Paleomagnetic  Remanence  and  Stratigraphy  of  Units 
--3-3.6  Ma.  Depicted  are  mean  CKRMs  for  each  lava  flow  sampled. 
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Figure  3.25  -  Paleomagnetic  Remanence  and  Stratigraphy  of  Renegade  Canyon  Area  Units. 
Depicted  are  mean  ChRMs  for  each  lava  flow  sampled.  For  localities  exhibiting  a  stratigraphy 
exhbiting  more  than  one  eruptive  unit,  stratigraphy  is  shown  to  indicate  eruptive  order. 
Stratigraphy  and  directions  are  colored  to  indicate  correlation  based  on  geochemistry  presented 
in  this  report.  Directiorvs  are  not  tilt-corrected. 
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Figure  3.26  -  Measures  of  Rotation  from  Lavas. 
Eruptive  unit  mean  directions  are  averaged  within 
locality  and  unit.  In  all  cases  except  F,  geographic 
roordinates  are  preferred  because  of  original  dips  of 
lava  flows  (e.g.  note  increased  scatter  in  A,  tilt- 
corrected).  (A)  Localities  within  Wild  Horse  Mesa  are 
not  rotated  relative  to  one  another.  (B&C)  Localities 
within  Wild  Horse  Mesa  are  rotated  relative  to  Cole 
Flat  area.  (D)  CHS  is  tilted  relative  to  Wild  Horse 
Mesa.  (E)  Results  suggesting  unaccounted-for,  large 
original  dips.  (F)  The  shallow  inclination  suggests 
that  secular  variation  is  not  fully  averaged.  (G&H)  No 
discemable  rotation.  No  tilt-correction  is  available  for 
these  latter  two  groups  of  localities, 
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■  Ar/Ar  presented  here,  with  2  a  error.  AM  other  geochronology  -  K/Ar  by  Duffield  et  al.,  1980.  Poiarity:  R  -  reversed,  N  -  normal. 


Figure  3.27  -  Volcano-stratigraphy  of  the  Study  Area.  Combined  use  of  stratigraphic  relations, 
geochronology,  and  paleomagnetic  polarity  reveals  this  sequence  of  major  eruptive  units. 
Some  minor  units  are  not  included.  Color  scheme  matches  that  used  in  other  figures. 
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Figure  4.1  -  Rotated  Blocks  of  Wild  Horse  Mesa.  This  cartoon  depicts  approxi¬ 
mate  boundaries  of  the  greater  Wild  Horse  Mesa;  the  Coso  Wash  fault  zone  and 
the  newly  named  Wild  Horse  Mesa  fault  zone.  Enclosed  are  blocks  oriented 
obliquely  to  both  boundaries  which  have  rotated  about  12°-20°  clockwise.  Three 
localities  for  which  we  have  direct  measures  of  rotation  are  shown.  Boundaries 
between  blocks  are  likely  sinistral-normal  faults.  Block  boundaries  curve  to  more 
north-south  orientations  towards  the  south  causing  a  space  problems  and  addi¬ 
tional  east  and  west  vergent  normal  faulting.  In  our  model  block  shapes  change 
to  accomodate  space  problems-  Likely  shape  change  is  indicated  at  block 
comers.  This  explains  the  observation  of  tihmst  faults  in  the  slickenline  data. 
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Analytical  methodsS  fresh  volcanic  rocks  were  chosen  for  ^Ar/^^Ar  geochronology 
HS-19,  Tdo,  RC-14,  RC-01  and  BC-18.  Rocks  were  crushed  and  sieved,  selected  fraction 
were  washed  in  water  and  in  an  ultrasonic  cleaner.  Using  mesh  40-60  (230  to  380  fim) 
fraction,  phenocrysts  of  olivine  were  removed  after  four  passes  through  a  magnetic 
separator  operating  from  0.5  to  up  to  1.7A.  The  presence  of  quartz  in  the  final  passes 
indicates  the  presence  of  xenocrysts  in  most  of  the  rocks,  which  have  basaltic 
composition  (e.g.  HS-19,  RC-01  and  BC-18).  In  order  to  avoid  the  spurious  effects  of 
possible  plagioclase  xenocrysts  we  have  separated  the  groundmass.  This  groundmass 
is  mainly  composed  of  fine  laths  of  plagioclase  (20-50  [im)  and  pyroxenes.  Transparent 
plagioclase  phenocrysts  from  the  same  samples  were  also  concentrated  by  magnetic 
separation.  Quartz  and  plagioclase  were  separated  from  the  non-magnetic  part  using 
heavy  liquid  (lithium  metatungstate  adjusted  at  density  =  2.645).  35  to  75  mg 
groundmass  and  plagioclase  samples  were  loaded  in  four  A1  disks  as  described  by 
Renne  et  al.  (1998)  and  irradiated  for  30  minutes  in  the  CLICIT  facility  of  the  TRIGA 
reactor  at  Oregon  state  University.  After  irradiation  25mg  of  plagioclase  and  90  grains 
of  groimdmass  for  each  sample  were  degassed  with  a  defocused  argon  laser  (Synrad) 
then  analyzed  with  a  MAP  215C  mass  spectrometer  at  the  Berkeley  Geochronology 
Center  using  methods  described  by  Renne  et  al.  (1998). 

Procediural  blanks  were  measured  between  eveiy  two  or  three  steps  and  typical 
values  were  (3.3  ±  0.5)x  10"^®  for  ^Ar,  (2.3  ±1.0)x  10'*®  for  ^Ai,  (5.3  ±1.0)x  10'^’  for  “Ar, 

(6.6  ±  1.0)x  10'^®  for  ®^Ar,  (2.26  ±  0.3)x  lO"^®  for  ®®Ar.  Neutron  fluence  (J)  was  monitored 
by  Fish  Canyon  sanidine  (PCs:  28.02  Ma,  Renne  et  al.  (1998))  grains  in  four  positions  in 
the  Al  disk.  The  J  value  for  each  disk  was  determined  from  individual  analyses  of  6 
separate  FCs  grains  in  each  position  fused  with  the  COj  laser  (total  of  24  grains).  The 
corresponding  J  values  calcifiated  with  the  total  decay  constant  of  Steiger  and  Jager 
(1977)  are  0.0001376  ±  0.0000007  for  Tdo  and  HS-19  plagioclases,  0.0001371  ±  0.0000010 
for  BC18,  RC-01  and  RC-14  plagioclases,  0.0001377  ±  0.00000006  for  Tdo  and  HS-19  for 
groimdmass  grams,  0.0001378  ±  0.00000007  for  BC18,  RC-01  and  RC-14  groundmass 
grains.  These  J  values  were  calculated  excluding  age  uncertainties  on  the  FCs  and  the 
decay  constant.  Mass  discrimination  was  monitored  by  automated  analysis  air  pipettes, 
interspersed  with  the  unknowns  was  1.0641  ±  0.0026  (2a  error)  per  mass  atomic  unit. 
Nucleogenic  production  ratios  from  K  and  Ca  were  (®®Ar/®^Ar)K=(1.22  ±  0.01)xl0r^ 
(®®Ar/®^Ar)ca=(2.19  ±  0.35)xl(l®.  Plateau  ages  are  defined  as  comprising  three  continuous 
steps  corresponding  to  a  least  50%  of  the  total  ®®Ar  released.  The  individual  fraction 
ages  should  agree  to  within  2a  of  the  'integrated'  age  of  file  plateau  segment. 


Results 

Groundmass:  Figures  1-5  show  the  ^Ar/®^Ar  incremental  heating  and 
ratio  spectra  and  inverse  isochron  plots  obtained  for  the  five  samples. 

HS-19:  The  spectrum  is  discordant  and  displays  a  dome  shape.  This  spectrum  is 
characterized  by  rising  ages  that  reach  old  apparent  ages.  The  last  35  %  of  the  spectrum 
correspond  to  lower  constant  ages.  Well  known  mechanisms  that  produce  discordant 
spectra  in  such  material  is  recoil  of  ®’Ar  during  irradiation  (Huneke  and  Smith,  1976; 
Forland  et  al.,  1993;  Sharp  et  al.,  1996).  In  the  summit  of  the  spectrum  (older  apparent 
age)  plagioclases  are  the  main  reservoir  of  ®®Ar.  As  a  consequence,  older  ages  are  found 
when  plagioclases  degas  unsupported  “Ar*  (Sha^  et  al.,  1996)  which  correspond  to 
intermediate  temperature.  The®’Ar  released  at  higher  temperatures,  principally  from 
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p)n:oxene  (higher  Ca/K  ratio,  see  Fig.  Bl),  without  corresponding  ^Ar*  result  in  low 
^parent  ages.  In  short  this  spectrum  shape  is  due  to  the  partial  redistribution  of  the 
from  the  plagioclase  to  pyroxenes  during  the  irradiation. 

Sample  Tdo  displays  a  saddle-shaped  spectrum  (Fig  B2).  Such  spectra  are 
characteristic  of  rocks  that  contain  excess  “Ar  (Dahymple  et  al.,  1975).  This  sample  has 
an  isotope  correlation  diagram  that  indicates  a  trapped  ^Ar/^Ari  ratio  of  379  ±  6.  This 
ratio  is  higher  than  that  of  atmospheric  at  the  95  %  confidence  level,  confirming  the 
presence  of  excess  argon.  This  excess  argon  is  generally  interpreted  in  lavas  flow  as  the 
result  of  incomplete  degassing  of  magmatic  Argon  (Sharp  et  al.,  1996).  We  interpret  the 
inverse  isochron  age,  2.83  ±  0.02  Ma,  as  the  maximum  age  for  the  Tdo  sample. 

RC-14  display  a  less  disturbed  spectrum  that,  after  a  slight  decrease  of  the  apparent 
ages  in  the  first  degassing  steps  displays  a  very  flat  spectrum  yielding  a  plateau  age  of 
2.62  ±  0.03  Ma  (Fig.  B3)  and  consistent  with  ^e  isochron  age,  2.607  ±  0.017  Ma.  The 
trapped  ^Ar/“Ar  ratio  (295.3  ±  1.1)  is  not  distinct  from  that  of  atmospheric  Ar. 

BC-18  and  RC-01  yield  ages  of  3.41  ±  0.04  (3.31  ±  0.03  isochron  age;  Fig.  B4)  and  1.37 
±  0.04  Ma  (1.30±  0.07  Ma  isochron  age;  Fig.  B5),  respectively.  Both  spectra  are  very  flat 
on  more  than  90  %  of  the  total  ®®Ar  released.  The  trapped  ^Ar/^Ar  ratios  are  not 
distinct  from  that  of  atmospheric  (297.3  ±  0.4  and  296.7  ±  1.1,  respectively). 

Plagioclase:  Four  samples  display  discordant  spectra  (HS-19,  RC-01,  BC-18,  and 
Tdo;  none  of  Ihese  are  shown).  The  complexity  of  these  spectra  could  be  attributed  to 
the  presence  of  plagioclase  xenocrysts  and  also  to  non-atmospheric  initial  (excess) 
argon,  which  contributes  to  the  apparent  high  ages  at  low  and  high  temperatures. 
Sample  RC-14  shows  a  plateau  age  according  to  our  criteria  defined  above  of  2.63  ±  0.04 
Ma  (97%).  This  age  is  not  distinct  from  that  one  obtained  from  the  coexisting 
groundmass  (compare  Figures  B6  and  B3). 
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Isochron  age  =  1 .30  ±  0.07  Ma 
40Ar/36Ar  Intercept  =  296.7  ±1.1 
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Isochron  age  =  3.25  ±  0,03  Ma 
40Ar/36Ar  Intercept  =  288.7  ±  1 .3 
MSWD  =  0.71 ,  Prob.  =  0.78,  n  =  18 
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Pliocene  volcanics 


□ 


Pre-Cenozolc  basement 


Figure  Cl  -  Lateral  Extent  of  Cole  Flat  Basalt  -  in  green.  This  unit,  comprised  of  at  least 
four  flows  near  its  source,  exhibits  normal  paleomagnetic  polarity  (declination  ~  335®,  incli¬ 
nation  ~  52®  ).  Question  marks  indicate  uncertainty  on  the  location  of  the  edge  of  the  unit. 
This  member  erupted  just  northeast  of  sampling  locality  CF  at  3.50  +/-  0.03  Ma.  See  main 
report  for  geochemical  composition  and  details  of  geochronology. 
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Figure  C2  -  Lateral  Extent  of  Black  Canyon  Flows  -  in  red.  This  unit,  comprised  of  at  least 
ten  flows  near  its  source,  exhibits  normal  paleomagnetic  polarity  (declination  ~  35®,  inclina¬ 
tion  ~  27®  ).  Question  marks  indicate  uncertainty  on  the  location  of  the  edge  of  the  unit.  This 
member  erupted  just  southwest  of  sampling  locality  BC  at  3.41  +/-  0.04  Ma.  See  main  report 
for  geochemical  composition  and  details  of  geochronology. 
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□  Pre-Cenozoic  basement 


Scale  Kilometers 


Figure  C3  -  Lateral  Extent  of  Haiwee  Spring  Basalt  -  in  aqua.  This  unit,  comprised  of 
nineteen  flows  at  locality  HS,  exhibits  normal  paleomagnetic  polarity  (declination  ~  320- 

330°,  inclination  ~  50-60°).  Question  marks  indicate  uncertainty  on  the  location  of  the  edge 
of  the  unit.  This  member  erupted  from  an  unknown  locality  north  of  the  sampled  localities 
indicated  around  3.32  +/-  0.1  Ma  (Duffield  and  Bacon,  1981).  See  main  report  for  geochemi¬ 
cal  composition. 
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Pliocene  volcanics 
Pre-Ceno2oic  basement 


Figure  C4  -  Lateral  Extent  of  Silver  Peak  Basalt  -  in  pink.  This  unit  exhibits  reversed 

paleomagnetic  polarity  (declination  ~  140-150°,  inclination  ~  -20° ).  Question  marks  indicate 
uncertainty  on  the  location  of  the  edge  of  the  unit.  This  unit  erupted  just  northwest  of  sam¬ 
pling  locality  SP  between  3.22  +/-  0.02  and  3.41  +/-  0.04  Ma  based  on  stratigraphic 
considerations.  See  main  report  for  geochemical  composition. 
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Pliocene  volcanics 
Pre-Cenozoic  basement 


Figure  C5  -  Lateral  Extent  of  Petroglyph  Canyon  Basalt  -  in  blue.  This  unit,  comprised  of 
at  least  six  flows  at  locality  LDC,  exhibits  reversed  paleomagnetic  polarity  (clockwise  rotated 

to  declination  ~  137°,  inclination  ~  -18°  ).  Question  marks  indicate  uncertainty  on  the  loca¬ 
tion  of  the  edge  of  the  unit.  This  member  perhaps  erupted  just  east  of  sampling  locality  FS  at 
3.22  +/-  0.02  Ma.  See  main  report  for  geochemical  composition  and  details  of 
geochronology. 
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Figure  C6  -  Lateral  Extent  of  QTbrj  Andesite  -  in  light  blue.  Correlation  of  this  unit  is 
tentative  because  of  equivocal  geochemical  correlation  and  the  lack  of  paleomagnetic  results  at 
locality  RCF.  Locality  RC  is  comprised  of  at  least  ten  flows  thought  to  be  QTbrj  and  exhibit¬ 
ing  an  unusual  reversed  paleomagnetic  direction  (declination  ~  250®,  inclination  ~  -75°). 
Question  marks  indicate  uncertainty  on  the  location  of  the  edge  of  the  unit.  This  member 
erupted  just  northeast  of  sampling  locality  RCF  at  2.62  +/-  0.03  Ma.  See  main  report  for 
geochemical  composition  and  details  of  geochronology. 
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Pliocene  volcanics 
Pre-Cenozoic  basement 


Figure  C7  -  Lateral  Extent  of  QTbrj  Andesite  -  in  orange.  This  unit  is  composed  of  a  maxi¬ 
mum  of  only  two  or  three  flows  (at  locality  RC).  QTbr2  exhibits  reversed  paleomagnetic 

polarity  but  variable  directions  (declination  ~  170-195°,  inclination  ~  -40  to  -50°).  Question 
marks  indicate  uncertainty  on  the  location  of  the  edge  of  the  unit.  This  member  erupted  near 
sampling  locality  RCE  at  1.37  +/-  0.04  Ma.  See  main  report  for  geochemical  composition  and 
details  of  geochronology. 
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Pliocene  volcanics 
Pre-Cenozoic  basement 


Figure  C8  -  Lateral  Extent  of  QTbr4  Basalt  -  in  plum  color.  Only  one  cooling  unit  of  QTbr4 
is  visible  at  the  localities  depicted.  QTbr4  exhibits  reversed  paleomagnetic  polarity 

(declination  ~  183°,  inclination  ~  -45°).  Question  marks  indicate  uncertainty  on  the  location  of 
Ae  edge  of  the  unit.  This  member  erupted  west  of  sampling  locality  RCD  at  an  unknown  time 
and  no  stratigraphic  constraints  are  available.  See  main  report  for  geochemical  composition. 
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